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t 

This  thesis  presents  a  procedural  sequence  for  processing  and 

assessing  the  integrity  of  experimental  strain  gage  data  taken  on  piping 
in  which  gage  rosettes  are  arranged  in  circumferential  patterns  at 
several  sections  of  a  piping  system.  The  method  of  least  squares  is 
used  to  obtain  "best1*  values  for  strain  per  unit  loading  for  each  gage 
element.  This  data  is  used  to  find  principal  stresses  and  values  of 
normal  and  shearing  stress  on  a  normal  cross  section  at  each  rosette 
location.  The  data  is  then  smoothed  to  obtain  normal  and  shearing  stress  as 
functions  of  angular  position  for  each  ring  of  rosettes.  A  computer 
program  has  been  devised  to  perform  these  operations  and  provide  a 
graphical  output  to  assess  degree  of  fit  with  experimental  data.  Finally, 
manual  computation  is  employed  to  study  the  equilibrium  of  various  seg- 
ments of  the  piping  in  evaluating  the  overall  integrity  of  the  data. 
This  procedure  and  program  have  been  used  in  evaluating  tests  performed 
on  the  main  seawater  piping  system  of  the  USS  Benjamin  Franklin  (SSBN 
640).  Some  results  of  the  dockside  hydrostatic  tests  are  presented 
herein  for  purposes  of  illustration. 
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1.  Introduction. 

ORIGIN  OF  IDEA  FOR  THESIS 

The  idea  for  this  thesis  was  conceived  while  the  writer  worked  at  the 
Electric  Boat  Division  of  General  Dynamics  Corporation  during  an  indus- 
trial tour  from  the  U.  S.  Naval  Postgraduate  School. 

While  working  on  the  analysis  of  strain  gage  data  taken  during  dock- 
side  hydrostatic  tests  and  sea  trials  of  the  USS  Casimir  Pulaski  (SSBN 
633 )»   it  became  evident  that  a  method  to  fully  analyze  the  main  seawater 
piping  system  from  the  strain  gage  data  did  not  exist.  The  only  tool 
available  at  the  time  was  a  computer  program  that  had  been  developed  by 
Electric  Boat  which  determined  the  two  principal  stresses  and  principal 
stress  axes  for  each  individual  set  of  rosette  readings. 

The  project  of  instrumenting  the  main  seawater  piping  system  on  the 
USS  Benjamin  Franklin  (SSBN  640)  and  taking  strain  gage  data  during  dock- 
side  hydrostatic  tests  and  sea  trials  had  been  assigned  to  Electric  Boat 
by  the  U.  S.  Navy  Bureau  of  Ships,  but  no  assignment  had  been  made  to 
analyze  the  data  taken.  The  writer  obtained  the  necessary  permission 
from  Electric  Boat  and  the  Office  of  the  U.  S.  Navy  Supervisor  of  Ship- 
building, Groton,  Connecticut,  to  use  the  data  in  connection  with  a  thesis 
at  the  U.  S.  Naval  Postgraduate  School  and  established  liaison  with 
personnel  at  Electric  Boat  connected  with  conducting  the  SSBN  6/fO  tests. 

Over  the  last  nine  months  Mr.  Carl  P.  Wilson,  Supervisor  of  the 
piping  section  of  the  Marine  Engineering  Design  Department  at  Electric 
Boat  and  his  main  assistant,  Mr.  William  R.  Hintz,  have  been  most  helpful 
and  prompt  in  forwarding  the  test  data  and  other  related  material  requested. 
Without  their  assistance  and  cooperation  this  thesis  could  not  have  been 
possible.  The  writer  would  like  to  express  his  sincere  appreciation  to 
Dr.  John  E.  Brock,  Professor  of  Mechanical  Engineering,  U.  S.  Naval 


Postgraduate  School,  for  his  continued  patience,  encouragement,  and  most 
capable  guidance  while  acting  as  faculty  advisor. 

The  intention  of  this  thesis  is  not  only  to  analyze  the  data  from  the 
SSBN  640, but  also  to  develop  a  procedural  sequence  for  processing  and 
assessing  the  integrity  of  strain  gage  data  that  could  be  utilized  in  the 
design  and  testing  of  future  submarine  seawater  piping  systems. 

,  HOW  DATA  WAS  OBTAINED 

Two  complete  sets  of  data  were  used  $&  developing  this  thesis. 
The  first  set  of  data  was  taken  on  7  August,  1965  during  the  dockside 
hydrostatic  test  of  the  main  seawater  piping  system  of  SSBN  640  in 
accordance  with  Electric  Boat  Division's  Propulsion  Plant  Test  Form 
Number  721-5  Revision  A,  The  second  set  of  data  was  taken  during  builder's 
sea  trials  on  19  and  20  September,  1965  in  accordance  with  Electric  Boat 
Division's  Sea  Trial  Form  Number  21.04  Revision  A.  However,  for  security 
reasons,  no  further  reference  is  made  in  this  theses  to  sea  trial  data. 

All  strain  gages  utilized  were  BLH  Type  FABR-25-12S9,  constantan  grid, 
phenolic  glass  base,  three  element  45°  superimposed  grid  rosettes*  More 
detailed  information  concerning  these  rosettes  may  be  obtained  from  BLH 
Specifications  Catalog  Bulletin  101  effective  September  1964. 

Gages  were  equally  spaced  around  a  cross  section  of  piping  at  30° 
from  each  other.  A  cross  section  witL  twelve  rosettee  will  be  referred  to 
as  a  ring  or  rosette  ring  throughout  the  remainder  of  this  thesis.  The 
rosettes  in  each  ring  were  oriented  as  shown  in  Figure  1,2*1.. ,  Ttapre  were 
eight  such  rings. 
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The  number  twelve  gage  is  located  facing  directly  up  for  all  rings 
except  ring  number  2,  in  which  it  faces  port.  The  approximate  locations 
of  the  rosette  rings  are  shown  in  figures  1.2.2,  1.2.3  and  1.2.4  which  are 
sketches  of  the  piping  system.  The  exact  location  of  each  ring  and  the 
piping  system  dimensions  are  not  given.  The  dimensions  may  be  found  on 
Bureau  of  Ships  Blueprints  3SBN  640  209  2117801  Rev.  H  and  SSBN  640  209 
2117802  Rev.  F  and  the  exact  ring  locations  may  be  obtained  from  either 
Electric  Boat  Division's  Propulsion  Plant  Test  Form  Number  721-5  Rev.  A 
or  Sea  Trial  Data  Form  Number  21.04  Rev.  A. 

Strain  gage  readings  were  taken  using  SR-4  Strain  Indicators  Type  N 
manufactured  by  Baldwin-Lima-Hamilton  Corporation. 

OBJECTIVES  OF  THESIS 

The  basic  purpose  of  this  thesis  is  to  develop  methods  of  assessing 
the  integrity  of  the  data  and,  incidentally,  to  reduce  the  data  to  a  form 
more  convenient  for  subsequent  analysis.  To  reach  this  goal  a  digital 
computer  program  has  been  developed  and  employed.  This  program  uses  as 
imputs  certain  constants  concerning  the  piping  material  and  sizes  and 
the  raw  strain  gage  data.  It  computes  the  two  principal  stresses,  the 
angle  measured  from  the  direction  of  flow  through  the  piping  of  the 
maximum  principal  stress,  and  the  shear  stress  at  each  cross  section  and 
stress  normal  to  each  cross  section.  In  addition,  smoothed  values  of 
normal  and  shear  stress  every.. 01  radian  around  the  circumference  of  the 
rosette  rings  are  calculated  and  plots  are  drawn  showing  both  smoothed  and 
actual  data  points.  Finally  the  program  computes  the  three  forces  and 
moments  at  each  rosette  ring. 

The  outputs  of  the  computer  program  are  analyzed  to  determine  the 
integrity  of  the  raw  data,  if  the  maximum  principal  stress  exceeds  yield 
strength,  the  agreement  between  smoothed  and  actual  data  values  of 
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normal  and  shear  stress  and  if  equilibrium  between  two  rosette  rings  is 
established. 

For  the  following  sections,  the  terms  experimental  data  and 
experimental  results  are  used  to  represent  the  actual  strain  gage  data 
taken  and  the  computed  results  from  this  data. 

Actual  dimensions  of  the  piping  system,  data,  and  actual  numerical 
results  have  been  eliminated  from  this  thesis  in  order  to  keep  the 
thesis  unclassified.  All  numerical  results  given  are  the  actual  results 
multiplied  by  arbitrary  constants.  However,  relative  values  of  the 
results  are  the  same  as  the  actual  results. 

NOMENCLATURE 


} 


a,b,c,d,     1    Constants  for  smoothed  value  of  normal  stress 


E  Young's  Modulus 

£i  Error 

e^  >  ^<r  Unit  vectors  in  the  ^    and  ^  directions 

F(0t),F,(0t)  Functions  of  the  angle  0 

Frk     •) 

FRr     }        Resultant  equilibrium  forces  in  the  x,y  and  z  directions 

Fr*     J 

Fx         Force  on  a  cross  section  of  the  main  seawater  piping  system 
in  the  direction  of  the  ship's  bow. 

Fy         Force  on  a  cross  section  of  the  main  seawater  piping  system 
in  the  vertical  direction  from  ship's  centerline. 

Fz         Force  on  a  cross  section  of  the  main  seawater  piping  system 
in  the  starboard  direction  from  ship's  centerline. 

f*'        7 

Fyt     V    Forces  in  the  x'  ,  y'  ,  and  z'  directions 
F*'    J 

i,j,k       Unit  vectors  in  the  x,  y  and  z  directions. 

k,l,m,n,  "1   Constants  for  smoothed  value  of  shearing  stress 

o,p,q   J 

xf7»^ix-'t7,tz/   Coordinate  systems  described  in  text 
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Mr* 

Mry     \        Resultant  equilibrium  moments  in  the  x,  y  and  z  directions. 

Mr* 

Mx         Moment  about  the  fore  and  aft  axis  of  the  ship  in  the 
direction  of  the  bow. 

My         Moment  about  the  vertical  axis  of  the  ship  in  the  upward 
direction. 

M?         Moment  about  the  port  and  starboard  axis  of  the  ship  in  the 
direction  of  starboard 

Mx'     ") 

My'     f        Moments  in  the  x'  ,  y'  ,  and  z'  directions. 

M*'  J 

r  Radius  of  the  piping. 

t  Pipe  wall  thickness. 

x  percent  of  test  depth. 

y  Strain  reading  at  any  depth 

[   3  Matrix 

[a]  a  7  x  7  matrix. 

[B]    ,      [b']]?x1  Matrices 
[6]  ,      [C«]  J 

€, ,  62  ,  63  Three  strains  of  a  three  element  90°  rosette. 

"^p  Maximum  principal  stress  angle 

p.  Poisson's  ratio. 

Z  Summation  sign 

<Tn\tm  Maximum  principal  stress 

(TmiK;  Minimum  principal  stress 

0"^  Stress  normal  to  a  cross  section  of  piping  in  the  direction 

of  flow. 

*r  Shear  stress  on  the  outside  surface  of  the  pipe  positive  in 

the  direction  of  increasing  rosette  numbers. 

X*%  Shear  stress  in  the  xz  plane 

tf  Angle  on  piping  cross  section  from  rosette  number  12. 

df  Stress  tensor  (dyadic). 
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2.  Theory1 

THEORY  FOR  STRAIGHT  LINE  FITTING 
It  is  assuned  that  the  strain  gage  readings  will  increase  linearly 
with  internal  pressure  which  in  turn  does  increase  linearly  with  an  in- 
crease in  depth.  This  assumption  is  based  on  the  following  facts: 

(a)  The  piping  material  exhibits  a  straight  line  stress-strain 
curve. 

(b)  The  system  is  not  loaded  above  the  elastic  limit  of  the 
material  as  verified  by  the  maximum  values  of  the  maximum  principal  stress 
at  every  rosette, 

(c)  The  deflections  are  small  as  verified  by  the  overall  change 
in  strain  gage  readings. 

(d)  It  is  assumed,  for  the  builder's  sea  trials,  that  hull 
deflection  inputs  to  the  piping  increase  linearly  with  depth.  For  dock- 
side  hydrostatic  tests,  there  was  no  such  input. 

The  other  main  effect  acting  on  the  system  is  the  weight  of  the  system. 
Since  the  system  is  always  filled  with  seawater  of  nearly  constant  density 
the  weight  effect  on  the  piping  system  will  remain  nearly  constant  and, 
therefore,  not  have  any  appreciable  effect  on  the  assumption  of  linearity. 
The  method  of  least  squares^  is  used  to  obtain  the  best  straight  line 
connecting  all  data  points  of  each  element  with  the  percentage  of  test 
depth  taken  as  the  independent  variable  and  the  strain  gage  readings 
taken  as  the  dependent  variable.  The  results  of  the  method  of  least 
squares  are  given  in  the  form  of  the  equation: 


"Greater  details  of  the  theoretical  relations  are  given  in,  or  may 
easily  be  inferred  from,  corresponding  subsection  of  section  3,  which  deals 
with  the  computer  program  used  to  perform  the  calculations. 

2 

Spiegel,  M.R.,  Schaum's  Outline  Series,  Theory  and  Problems  of 

Statistics  (New  York;  Schaum. Publishing  Co.,  1961)  p. 220 
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Y  -  a  +  hx  2.1.1 

where  a  is  the  zero  percent  test  depth  intercept  and  b  is  the  slope  of 
the  best  straight  line  connecting  all  data  points  from  the  surface  to 
test  depth.  The  slope,  b,  is  also  proportional  to  change  in  strain  per 
percent  change  in  test  depth.  The  slopes  are  the  values  of  the  indi- 
vidual element  strains,  6^,  (%>»  an(*  ®3>  used  in  the  subsequent  develop- 
ment of  the  theory. 

THEORY  FOR  SOLVING  PRINCIPAL  STRESSES  AND  ANGLES 
The  equations  defining  the  principal  stresses  in  terms  of  the 
strains  measured  by  a  rectangular  rosette  are3j 


2 


fLtfiL+     i"  'VC6.-602-   +r2€2-(6,f6^fx?    2.2.1 
I  l-A        1+/JL  J 

CTmiM    *i.fi^a.-  J_VC6,^3)l+L"26l-(€(+f3)]1   ]     Z2-2 
21  l-XJL        i+u  J 


where  the  elements  of  the  rosette  are  numbered  as  shown  in  figure  2.2.1 

3 

a^  5 —  Direction  of  algebraically 
greater  stress. 

,eP 


Figure  2.2.1 
The  direction  of  the  maximum  principal  stress  from  element  number  one 
of  the  rosette  in  terms  of  the  measured  strains  iss 


Op  =  JL  <xrc  +an      26z  -  (6. +  6*) 
2  e,-e* 


2.2.3 


3 

-'The  derivation  of  equations  2.2.1,   2.2.2,  and  2.2.3  is  given 

in  Perry  and  Lissner,   The  Strain  Gage  Primer,   pages  136-138 
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THEORY  FOR  EXPERIMENTAL  VALUES  OF  CT*  #  t*r 
Looking  down  on  the  surface  of  the  piping,  the  stress  picture  looks 
like  figure  2.3 .1. 


direction  of  flow 


Figure  2.3.1 
It  is  known  that  0^  which  is  perpendicular  to  the  plane  of  (X  and  Gl 
is  equal  to  zero. 

To  determine  the  stresses  on  a  plane  perpendicular  to  the  x-y  plane 
whose  trace  is  mm,  tensor  theory  is  used.  Temporarily  use  coordinates 
^"P^and  z  to  determine  CT/^  t   Lxy  andlTxe.  (The  direction  z  is  per- 
pendicular to  the  plane  shown  in  figure  2.3.1  with  unit  vector  k  .) 
From  tensor  theory,  the  stress  tensor  (dyadic)  may  be  written: 

^=  Gm/tx  CfCf   +  (Tm.w  e-»  e 7  2.3.1 

<Tm=  T-^-I  2 .3.2 

txy-Tv^.J  2.3.3 

tx*  =  Iv«f-£  2.3.4- 

Since  the  term  f-^   appears  in  equations  2.3.2,  2.3.3,  and  2.3.4,   solve 
for  it  first. 

=    CknAx  Cc?S  Op  ee  -  CTmiM  Sinj  B-p  e ^  2-3-5 

Substituting  equation  2.3.5  back  into  equations  2.3.2,  2,3.3,  and  2.3.4 
results  in  the  following  equations: 

(Tm=    Ornftx   COS  Op  (Cc'l)    '  Cm.io    SlWOp   (<?»•  7) 

=  OmRx  C0S2O-p  4  OmiM    Siw2-Op  2.3.6 
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=   Om*x  COSOp  SliOBp  -   CTm in    S»N)©p   COSOp 
=     /^((Tmax  -  (Tmiw)  SJN  2  Op  2.3-7 

txa  =  (Tmftx    CosG-p(ee-  k)  -  dmi*   SlfJ-Op  (  e^>  •  k) 
=  (flm**  COSa-p)  O  -(Om»N  SiM-Gp)  O 
=  O  2.3.8 


TKE0&Y  FOR  SMOOTHED  VAXUSS  OP  NOffifAL  AMD  SHEARING  STRESS 
In  developing  the  theory  for  the  smoothed  values  of  On  and  txy 
it  has  been  assumed  that  the  radius  of  the  piping  has  remained 
constant  throughout  the  entire  piping  system.     This  assumption  is  based 
on  the  diameter  readings  taken  by  GD/EB  of  rosette  rings  two  through 
eight.     These  readings  indicated  that  the  diameter  of  the  piping  around 
a  rosette  ring  remained  constant  enough  to  use  a  constant  value  of  radius 
for  calculations  in  this  thesis. 

Since  a  previously  developed  theory  was  not  available,  it  seemed 
reasonable  to  assume  that  On  and  t*Y  vary  only  with  0  in  the  following 
Fourier  Series: 
Cn-  Gl+bs!N/<#-f-Ccos<£-h  d.S/N/20  +•  e  COS  20  +  £  Sin  30 

+  cjcos30+- =   F(0)  2. 4-1 

txy  -    k+   ISIN0-I-    m  COS  <f>  +■  71  Si"  20  -h    O  COS  Z<j/>  +■  p  S/AJ  30 

+  ycos  30+">>  -  F, (0)  2.4-2 

One  of  the  decisions  that  had  to  be  made  concerned  the  number  of 
terms  of  the  Fourier  Series  to  use.  Since  curves  showing  the  exact 
shape  of  the  stress  distribution  around  the  circumference  did  not  exist, 
there  was  no  criterion  available  to  which  curves  developed  could  be 
compared.  To  assist  in  deciding  how  many  terms  to  use,  the  solution  was 
investigated  using  five,  seven,  nine,  and  eleven  terms.  A  plot  of  the 
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experimental  values  and  the  smoothed  values  of  On  for  five,  seven,  nine, 
and  eleven  terms  for  rosette  ring  number  three  is  shown  in  figure  2.4.1. 
A  decision  was  made  to  use  seven  terms  for  the  following  reasons: 

(a)  A  certain  amount  of  redundancy  was  necessary  to  permit 
good  least  squares  evaluation. 

(b)If  higher  coefficients  were  really  of  significance,  more 
than  twelve  rosettes  would  be  required  at  each  rosette  ring  to  obtain 
more  than  seven  coefficients  whose  accuracy  could  not  be  doubted. 

(c)  For  equilibrium,  only  the  first  three  coefficients  are  used. 
Therefore,  it  is  necessary  to  have  these  evaluated  as  well  as  possible. 

The  truncated  forms  of  equations  2.4.1  and  2.4.2  shown  below  will  be 
used  in  the  subsequent  development  of  theory. 
On*  a  +   bs»W(£  +  CQOS(/>  +  d.S»W  20  +-  6  COS  2^  f  fs/w3^ 

-i-cjcos3<#=  F(0)  2.4.1a. 

t<Y  -  ^  +•    Is  IW0   +-7UCOS0  +TlS(w20f   O  COS  20  +■  pS/W30 

4-cjcos30  s  P.C0)  2.4.2a. 

The  smoothed  values  of  0*will  be  discussed  before  proceeding  with  the 
smoothed  values  of  t*y. 

There  are  N=12  stations  around  the  circumference  of  the  pipe  for 
which  there  are  experimental  values  of  0"^  ,  i.e.,  there  exists  a  (Tc 
corresponding  to  F(0[  ),  i^ij—^N.  0~t  is  considered  to  be  the  experi- 
mental value  if  F( 0j)  is  considered  to  be  the  smoothed  value.  It  is 
assumed  that: 

error  =  6c   =  &C  -  F  (0c)  2.4-3 

It  is  necessary  to  minimize  the  sum  of  the  square  of  the  errors  in  order 
to  have  smoothed  results  and  experimental  results  within  reasonable 
agreement  ,j  Therefore, 
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FIGURE  2.4.1 
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2.4.4 


is  to  be  minimized.  This  means: 

dO- 

$-° 

if    ^ 

"Si"  7 

Equations  2.4.4  are  solved  in  the  following  manner: 

ec=(T:-  F(tfi) 

=  CTc-[a+bsiN0  +  ccos<£-i-clsiM2^>4ecos296  +- £sin30+c\cos30]  24.5 

^£=  -ZJL  (n  ^F«&)  +  2F<*t)AEfflSd-  O  2.46 

where  K  =  a,  b,   c,  d,  e, ,f,  g 

By  substituting  the  various  values  of  K  in  equation  2.4.6  the  following 
equations  are  obtained: 

3a 

_ag  ^   -  25  gv  siKJC^t'  +■  21  F(0c)  sin0,:=  O 
£b 

££^  -  21CU-SIM  20c  +  2<ZF(&)sin20c'  =  O 
|£=  -  2  £07  cos  20;  -h  2£FW>c)cos2<j6;  =  0 


2.4.T 
2.4- S 
2.4.9 
2.4. IO 

2.4.H 
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2£»-  22<ncos3<^t-+  22t((t>c)C0S3<t>r=O  2.4.13 


Shifting  the  first  terms  to  the  right  hand  side  and  substituting  equation  2-4-Icl 
for  F(<$t)  in  equations  2.4.7  through  2.4.13,  the  results  displayed  in 
matrix  form  are? 

where: 


M 


M- 


[o]= 


■ 

Ss0t 

is2$ 
2c20i 

£s30c 
*c30i 


SC0; 

.  riff 

0t-  Is20c0; 
0L    £c20p0i 

ft    £c30(c0c 


Ia20;  2c20t 

Is0s20t  £s0p20; 

2c#s2ft  lcftc20; 

£s22fc  Ss20c2# 

£c20fe20c  2c2^t 

2s3jte20c  2s30c20;      ls230-         Is30|c3ft 

ic30s2ft  Sc3ftc2ft      Ic30s30-     Sc230i 


Zs30f  Sc30t 

2s{fe30L  ^s0c30c 

£c^30(  *Zcfe30c 

2.s20s30c  ^s20p30,: 

2c21te3#  Sc2ftc30i- 


2.4.15 


(In  this  matrix  only,  s  is  used  for  sin,   c  is  used  for  cos.) 


a 

b 
c 
d 
e 
f 
g 

£GTisin0i 

lfl"i.cos0t 

S(Flsin20j 

2d[cos20r 

£<T  |  sin30£ 

l(Ticos30i. 


2.4.16 


2.4.17 
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When  the  number  of  stations  around  the  circumference  is  even  and  the 
stations  are  equally  spaced,  it  can  readily  be  shown  that  all  non-diagonal 
elements  of  the  A  matrix  in  equation  2.4.15  become  zero.  In  this  case  the 
seven  equations  vdth  seven  common  unknowns  become  seven  equations,  each 
with  only  one  unknown,  as  follows: 

a  A/  =£a~<k 

C  £ Cosx(j>  i  =  5  (TV  COS 1 0C 

ci$siwl2<tfi*  -Savs'rw*^  \    2.4.13 

e  £  Cos7  2(/;   =  ^GV  cos  v  20t* 
££  SIN)2-  30i  =   &TC  S/aj  l  30<' 

cj  Ecos1  3<A*  -  5  or  <°os z  30 c" 

By  a  similar  derivation,  the  following  equations  are  established  for 
the  seven  coefficients  used  in  determining  the  smoothed  values  of  Uy. 

[A]  [b']  -  [C']  2.4.19 

where  [a]  is  the  same  as  shown  in  equation  2.4.15* 


M 


M- 


k 

1 
m 
n 
o 
P 

q 


lt|sin0| 

5t-'  cos0t 
Itj  sin20r 
£'t;  cos20c 
^T •:.  sin30: 
1't  c  cos30L 


2.4.20 


2.4.21 


Equations  2.4.17  become: 


kN  =  itj 

l£sin0  =  ?t,sin^L 

mlcosp  =  5iicos0j 

n2sin220  ■  £tisin20^ 

o2cos220  =  f>ticos*20j 

p£sin230  =  i-tisin230; 

q£  cos23^  -  Sticos^jfe 
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2.4.22 


Finding  the  seven  unknowns  for  On  and  the  seven  unknowns  for  Txyby  use  of 
equations  2.4.18  and  2.4.22  will  be  referred  to  as  the  short  method  in 
the  remainder  of  this  thesis.  Finding  the  unknowns  for  0"n  and  tx-yby  use 
of  equations  2.4.14  and  2.4.19  will  be  referred  to  as  the  long  method  in 
the  remainder  of  this  thesis. 

Once  the  two  sets  of  seven  unknowns  have  been  found  by  either  the 
short  or  long  method,  they  are  substituted  into  equations  2.4.1a  and 
2.4.2a  to  determine  the  smoothed  values  of  On  and  T*y  . 

THEORY  FOR  EQUILIBRIUM  CALCULATIONS 

Equilibrium  of  the  system  was  established  in  several  steps  .   The 
first  step  was  to  calculate  the  sum  of  the  forces  and  moments  at  each 
rosette  ring  cross  section.  The  x'  direction  was  taken  as  the  direction  of 
flow  in  the  pipe,  the  j*   direction  as  the  direction  from  the  center  of 
the  pipe  toward  rosette  number  12,  and  the  zf   direction  as  the  direction 
from  the  center  of  the  pipe  toward  rosette  number  3>  as  shown  in  figure 
2.5.1. 


*-x 


Figure  2.5.1 
The  equations  to  obtain  the  sum  of  the  forces  and  moments  are: 
Fx'  =-tl  rfod.0  *  rtj    (a+bsiN0  ^-ccos0+  dsiw  20-i-"  . )  d0 
=  2TTart  °  2.5.1 

Fy*  -t  \     rSiU^Txy  d.0  =  -rt  \  StM0(k+\siM&  +  rc\cos0  +  >>  <<)d$ 


=  -TTlr-t 


2.5.2 
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*/7  .    ,  ,     ^2/7- 


=    77/vr-f  2.5.3 

=    ZTTkr-t  2.5.4- 

•/a  Jo 

e-TrLr^t  2.5-5 

/a  J  & 

=  -  rrcrijr  2.5.6 


The  second  step  in  the  equilibrium  calculations  is  to  reorient  each 
individual  cross  section  coordinate  system  to  a  standard  coordinate 
system.  The  standard  coordinate  system  aligns  the  x  direction  with  the 
fore  to  aft  axis  of  the  ship,  the  y  direction  with  the  vertical  upward 
axis  of  the  shipj  and  the  z  direction  with  the  port  to  starboard  axis  of 
the  ship.  The  reorientation  is  obtained  using  actual  dimensions  and 
angles  taken  from  the  system  blueprints  and  standard  coordinate  re- 
orientation methods. 

The  third  step  in  the  equilibrium  establishment  consists  of  making 
the  segments  between  rosette  rings  into  free  bodies  and  summing  the 
forces  and  moments  acting  on  the  segment.  In  addition  to  the  forces  and 
moments  calculated  in  equations  2.5.1  through  2.5.6,  the  forces  due  to 
the  fluid  in  the  pipe  segment  are  added  to  obtain  the  equilibrium  forces 
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and  moments.  The  signs  of  the  forces  and  moments  at  the  end  of  the 
segment  in  which  the  fluid  enters  are  reversed  because  the  values  cal- 
culated in  equations  2.5.1  through  2.5.6  are  the  forces  of  the  adjacent 
segment  acting  in  that  segment  and  not  the  forces  of  the  incoming  segment 
acting  in  the  adjacent  segment  as  are  desired  for  equilibrium  calculations. 
Theoretically  the  sums  of  the  forces  and  moments  on  each  individual 
segment  should  be  zero. 

Only  step  one  utilizes  the  computer  program  developed.  Since  each 
rosette  ring  is  oriented  differently  with  respect  to  the  standard 
coordinate  system  chosen,  more  time  would  be  required  to  adjust  the  program 
to  calculate  the  different  reorientations  for  equilibrium  calculations 
than  to  calculate  the  reorientations  by  hand.  The  same  applies  to  step 
three.  Therefore,  steps  two  and  three  are  performed  by  hand. 
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3.  Explanation  of  Computer  Program 
$gi      i  .   INTRODUCTORY  COMMENTS 

This  section  presents  a  full  explanation  of  the  digital  computer 
program  developed  for  evaluating  the  strain  gage  data.  A  listing  of  the 
various  symbols  used  in  the  program,  but  not  defined  in  the  program 
explanation,  is  found  in  Appendix  B,  A  complete  listing  of  the  program 
appears  in  Appendix  A.  The  statement  numbers  referred  to  in  this  section 
are  the  six  digit  letters  following  each  line  of  the  program  listing,, 

The  program  language  employed  in  this  program  is  CDC  FORTRAN  60. 
All  computer  runs  were  made  on  the  Control  Data  Corporation  1604  computer 
and  associated  equipment  located  at  the  U.  S.  Naval  Postgraduate  School 
at  Monterey,  California. 

The  program  developed  is  divided  into  12  parts.  The  following  sub- 
sections describe  these  parts. 

PART  1,  DATA  READ  IN 

Part  1  of  the  program  reads  in  various  quantities  that  remain  constant 
during  the  entire  program.  These  quantities  most  likely  would  vary  if 
used  for  another  piping  system,  but  would  still  remain  constant  for  that 
particular  system. 

The  quantities  read  in  are  N,  the  number  of  experimental  depths  at 
which  readings  were  taken;  E,  Young's  modulus  for  the  pipe  material  in 
millions  of  pounds  per  square  inch;  n    Poisons  Ratio  for  the  pipe 
material  used;  TH(l)  through  TH(8),  the  average  wall  thicknesses  for  rings 
one  through  eight;  RA(l)  through  RA(8),  the  pipe  radius  for  ring  one 
through  ring  eighth;  and  Z,  the  angle  between  two  strain  gage  rosettes  of 
each  ring  in  radians.  Also  included  in  this  part  are  the  commands  which 


^These  values  of  radius  and  thickness  were  obtained  by  averaging 
the  thicknesses  and  radii  at  the  various  measurement  points  for  each  ring. 
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initialize  all  AM,  BM,  and  CM  registers  at  the  start  of  calculations 
for  each  ring. 

Statement  number  000280  indicates,  the  number  of  rings  of  strain  gage 
rosettes  for  which  data  is  to  be  computed.  This  is  indicated  by  the 
number  of  times  the  nD0  LOOP"  is  used. 

PART  2,  COMPUTATION  OF  BEST  SLOPE 

Part  2  of  the  program  reads  '/he  actual  strain  gage  readings  in  micro- 
inches  per  inch  from  data  cards  at  the  end  of  the  program  and  by  the 
method  of  least  squares  determines  the  slope  of  the  line  connecting  the 
data  points  from  0$  to  100$  of  test  depth.  The  slope  is  given  by  the 
value  of  B(M)  for  each  gage.  The  actual  strain  at  zero  percent  depth  is 
given  by  the  calculated  intercept,  identified  in  the  program  by  AZ.  The 
zero  percent  depth  intercept  is  not  used  in  this  program,  however. 

An  indication  of  the  accuracy  of  the  data  can  be  determined  in  this 
part  of  the  program.  One  of  the  quantities  calculated  in  this  part, 
identified  in  the  program  by  DIFFS(I),  gives  the  difference  between  each 
actual  strain  gage  reading  and  the  calculated  best  line  fitting  value  at 
each  reading  depth.  The  magnitude  of  DIFFS(I)  will  indicate  the  accuracy 
of  the  data. 

In  (future  applications  DIFFS(I)  could  be  used  to  expand  the  program 
to  accept  or  reject  each  and  every  data  point  by  inspecting  the  magnitude 
of  DIFFS(I)  compared  to  the  values  of  DIFFS(I)  for  a  particular  gage  at 
other  depths.  The  program  at  present  does  not  reject  any  individual  data 
points.  If  one  or  more  points  were  to  be  rejected,  the  program  could 
be  altered  to  compute  a  modified  slope  and  intercept. 

The  slope,  B(M),  in  this  program  is  given  in  strain  per  percent  of 
test  depth. 

Part  2  of  the  program  is  run  three  times  before  going  on  to  part  3  of 

26 


the  program  in  order  to  nave  all  three  strains  available  at  once  when 
computing  the  principal  stresses  for  each  rosette. 

PART  3,  COMPUTATION  OF  PRINCIPAL  STRESSES  AND  AXES 

Part  three  of  Lhe  program  utilizes  the  strain  per  percent  of  test 
depth  calculated  in  part  two  of  the  program  for  each  element  of  the 
rosettes  to  determine  the  two  principal  stresses  and  direction  of 
maximum  principal  stress  taking  the  direction  of  flow  within  the  pipe  as 
the  zero  angle.  The  values  of  the  two  principal  stresses  are  given  in 
pounds  per  square  inch  per  percent  of  test  depth. 

Statements  000860,  000870,  and  000880  check  to  see  if  any  of  the 
individual  element  strains  are  zero.  If  any  individual  strain  is  zero, 
the  program  jumps  to  statement  001020  which  sets  the  two  principal 
stresses  and  angle  equal  to  zero. 

Statements  000890  through  001000  calculate  two  principal  stresses 
using  the  formulas 


2U-m  L+M-  J 


Tli-M-   1+M-  J 

and  the  angle  using  the  formula 

An^Ic  =_L  arc  law  2ei-  ((^  +  €3) 

Part  2  and  part  3  of  the  program  up  to  statement  001050  is  run  twelve 
times  before  continuing  further,.  The  number  of  times  these  two  parts  are 
run  is  determined  by  the  number  of  rosettes  per  rosette  ring.  For  the 
data  evaluated,  there  were  twelve  rosettes  per  ring. 

Statements  001060  through  001090  orient  the  angle  of  the  maximum 
principal  stress  with  respect  to  the  direction  of  flow  through  the  pipe. 
This  is  only  necessary  when  one  of  the  outside  legs  of  the  rosette  is  not 
oriented  in  the  direction  of  flow.  The  rosette  orientation  as  indicated 
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in  figure  1.2.1  has  caused  the  necessity  in  this  particular  case  for 

0 

including  statements  001060  through  001090. 

A  print  out  of  values  of  O^m^  ,Grnnsj  and  -Op  for  all  twelve  rosettes 
in  a  ring  is  called  for  in  statements  001100  and  001110. 

PART  4>  DECISION  REGARDING  MJSFHOD  OF  SOLUTION 

Part  four  of  the  program  determines  whether  the  short  or  long  solution 
to  solve  for  the  smoothed  values  of  On  and  Xx.y  must  be  used. 

In  subsequent  sections  of  this  thesis  the  use  of  the  word  inoperative 
may  have  one  of  two  meanings.  In  one  meaning,  it  is  used  to  represent  any 
rosette  element  for  which  data  was  not*  recorded-*.  In  the  other  meaning, 
it  is  used  to  indicate  any  rosette  element  with  a  zero  slope  as  determined 
by  the  method  of  least  squares.  The  second  meaning  was  a  decision  made 
by  the  writer  based  on  a  close  evaluation  of  the  data  without  its  use. 

When  any  element  of  a  rosette  is  inoperative,  the  value  of  CTrnAx.  , 
calculated  in  part  three,  for  that  rosette  is  set  equal  to  zero.  This 
part  of  the  program  checks  to  see  if  any  value  of  Omfixfor  a  ring  is  zero. 
If  any  one  value  of  (TWiax  is  zero,  the  program  jumps  to  the  long  method  of 
solution.  If  none  of  the  values  of  Cmn*  for  a  ring  is  zero,  the  program 
goes  to  the  short  method  of  solution. 

PART  5,  DETERMINATION  OF  FOURIER  COEFFICIENTS  -  SHORT  METHOD 
Part  five  of*  the  program  determines  the  values  of  the  coefficients 
used  in  the  determination  of  the  smoothed  values  of  G"n  and  0TxY  by  means 
of  the  short  method  as  defined  in  section  2  of  this  thesis.  The  computer 
space  names  are  used  in  both  parts  five  and  six  in  determining  the  same 
items. 

5When  data  has  not  been  recorded  for  a  rosette  element t     the 
strain  readings  entered  on  the  program  data  card  for  that  element  are 
zero.  This  is  done  in  order  to  keep  proper  data  card  sequence  in  the 
program, 
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Statement  001230  changes  the  angle  ^from  degrees  to  radians  while 
statements  001240  through  001260  solve  for  the  experimental  values  of  GVj 
and  txy  •  Phi,  statement  number  001280>  is  the  angle  of  the  rosette  from 
rosette  number  12.  Statements  001290  through  001400  compute  the  values 
of  the  sines,  cosines,  sines  squared,  and  cosines  squared  of  the  single 
angle  phi,  twice  the  angle  phi  and  three  times  the  angle  phi.  The  value 
of  AM(1)  in  statement  001410  indicates  the  number  of  rosettes  in  each 
ring  while  statements  001420  through  001470  sum  up  the  squares  of  the 
sines  and  cosines  of  the  single,  double  and  triple  angles,  phi.  These 
last  seven  statements  are  the  main  diagonal  elements  of  the  7x7  matrix 
shown  in  theory  (equation  2.4.15).  Statements  001480  through  001610 
calculate  the  values  of  the  two  7x1  matrices  shown  in  theory  (equations 
2.4.17  &  2.4.21).  Statements  001620  through  OOI65O  change  the  values  of 
Gn  and  X*.y   into  values  needed  in  making  the  polar  plots.  The  last  number 
of  the  "DO  LOOP"  which  connects  statements  001220  through  001660  is 
determined  by  the  number  of  rosettes  in  the  ring. 

Statements  001690  through  001780  are  name  changing  statements. 

Statements  001790  through  001850  solve  the  seven  equations  for  the 
seven  unknown  coefficients  for  the  smoothed  values  of  On  while  statements 
001860  through  001920  solve  the  seven  equations  for  the  seven  unknown 
coefficients  for  the  smoothed  values  of  t*Y  .  These  14  equations  are 
equations  2.4.18  and  2.4.22  in  section  two  covering  theory.. 

The  remaining  statements  in  this  part  are  all  print  and  format  state- 
ments to  give  a  print  out  of  various  items  computed  in  this  part. 

When  part  five  is  used  in  the  program,  statement  002010  is  used  to 
bypass  part  six  of  the  program  and  go  directly  from  part  five  to  part 
seven. 
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PART  6,  DETERMINATION  OF  FOURIER  COEFFICIENTS  -  LONG  METHOD 

Part  six  of  the  program  determines  the  values  of  the  coefficients  used 
in  the  determination  of  the  smoothed  values  of  On  and  X*.^   by  means  of  the 
long  method  as  defined  in  section  2  of  this  thesis. 

Statement  002080  is  used  to  determine  if  any  rosette  in  the  ring  is 
inoperative.  If  any  rosette  is  inoperative,  statements  002090  through 
002150  are  used  and  the  program  goes  to  the  end  of  the  "DO  LOOP".  This 
is  used  to  prevent  plotting  an  experimental  value  of  On  and  Xxy  on  the 
output  plots  when  the  rosette  is  inoperative. 

Statements  002160  through  002290  are  identical  to  statements  001230 
through  001360  of  part  five.  Statements  002300  through  002480  compute 
the  combinations  of  sines  and  cosines  of  the  single,  double  and  triple 
angle  phi.  The  items  identified  by  AM(l)  through  AM(28),  statements 
002490  through  002760,  are  the  various  elements  of  the  7x7  matrix  shown 
in  section  two  (equation  2.4.15).  Statements  002770  through  002830  are 
the  elements  of  the  7x1  matrix  shown  in  theory  (equation  2.4.17)  while 
statements  002840  through  002900  are  the  elements  of  the  7x1  matrix 
shown  in  theory  by  eolation  2.4.21.  Statements  002910  through  002940 
are  identical  to  statements  001620  through  001650  in  part  five  of  the 
program.  The  last  number  of  the  "DO  LOOP"  connecting  statements  002070 
through  002950  is  once  more  determined  by  the  number  of  rosettes  in  each 
ring. 

Statements  002980  through  003460  identify  all  elements  of  the  7x7 
matrix  individually  (Equation  2.4.15)  as  can  be  noted  by  the  double 
subscripts.  In  turn,  statements  003500  through  003520  do  the  identical 
identification  for  the  two  7x1  matrices  (Equations  2.4.17  and  2.4.21). 
Statements  003530,  003540  and  003550  are  the  first  three  items  needed  in 
the  CALL  GAUSS  2  statements.  The  two  CALL  GAUSS  m statements,  numbers 

r 
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003560  and  003610,  call  the  subroutine  GAUSS  2  into  use.  GAUSS  2 
solves  for  the  unknowns  in  the  matrix  equations  and  is  used  twice  to 
solve  for  the  unknown  coefficients  of  both  CTw  and  Uy« 

The  various  other  statements  in  this  part  of  the  program  are  print 
and  format  statements  which  print  out  various  quantities  computed  in  this 
part. 

PART  7,  PREPARATION  FOR  GRAPHICAL  OUTPUT 

Part  seven  of  the  program  is  directly  related  to  the  subroutine  DRAW. 
It  is  merely  a  series  of  statements  which  set  up  the  labels  placed  on 
figures  drawn  by  the  program.  More  information  on  the  titles  and  the 
requirement  for  them  can  be  obtained  by  consulting  Appendix  D.  State- 
ments 003750  through  003900  are  a  computed  GO  TO  series  for  numbering  the 
rings  appropriately  as  computer  draws  the  figures. 

PART  8,  COMPUTATION  OF  SMOOTHED  VALUES  OF  AXIAL  AND  SHEARING  STRESS 

Part  eight  of  the  program  computes  the  smoothed  values  of  Gw  and  tx.Y 
around  the  entire  pipe  diameter  and  converts  the  values  into  polar 
coordinates  for  plotting  purposes.  The  360  degrees  of  the  circle  of  the 
pipe  circumference  is  divided  into  628  segments  for  plotting  purposes. 
Since  the  computer  plotter  used  for  figure  drawing  can  only  draw  vertical, 
horizontal  and  forty-five  degree  lines,  it  was  necessary  to  divide  the 
pipe  circumference  into  as  many  convenient  increments  as  possible.  The 
value  628  was  chosen  since  the  computer  utilizes  radians,  not  degrees  for 
angle  s;  and  628  is  a  convenient  number  with  which  to  divide  27T  radians 
and  900  is  the  limiting  number  of  points  that  can  be  plotted  by  the  sub- 
routine DRAW. 

Statement  003970  divides  2TT  by  628.  Statement  004000  computes  the 

angle.  Statements  004010  through  004060  compute  the  sines  and  cosines 

of  the  angles  needed  in  the  formulas  f or  On  and  /hcY  •  Statements  004070 

*  ■ 
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and  004080  compute  the  smoothed  values  of  0\i  using  the  formula 
On=GL+6s/m0  +  ccos<j6  +  cIs</v20  +  ecos20  -t-^SiN  30  +  Q  Cos  30  • 

Statements  004090  and  004100  compute  the  smoothed  value  of  txy using  the 
formula 

txY=  k  +  Istv0  +  mcos<f>+n$iuz0  +  a  cos  20  -t  p  s/aj  30  +  ycos  30 . 
Statements  004110  through  0041&0  transform  the  values  of  CfNand  txy  into 
values  that  can  be  utilized  in  polar  plotting.     Statements  004110,   004120, 
004150,  and  004160  are  used  for  drawing  reference  circles  on  the  figures. 
These  reference  circles  indicate  the  zero  stress  ring.     Any  reading  inside 
the   zero  stress  circle  for  Ok/ can  be  taken  as  a  compression,  while  any- 
reading  outside  the  zero  stress  circle  for  On  can  be  taken  as  a  tension. 
For  the  plots  of  %^,&nj  reading  inside  the  zero  stress  circle  can  be 
taken  as  a  negative  shear,  while  any  reading  outside  the  zero  stress  circle 
can  be  taken  as  a  positive  shear.     A  positive  shear  is  defined  as  a  shear 
caused  by  the  portion  downstream  of  the  rosette  ring  acting  in  the  direction 
of  increasing  rosette  numbers  on  the  upstream  portion.     A  negative  shear 
is  defined  as  a  shear  caused  by  the  portion  downstream  of  a  rosette  ring 
acting  in  the  direction  of  decreasing  rosette  numbers  on  the  upstream 
portion.     The  number  100,   chosen  for  plotting  the  zero  stress  circle  for 
(Jn  and  40,   chosen  for  plotting  the  zero  stress  circle  for  fcy,  were  chosen 
to  keep  all  plotted  values  positive  and  to  show  as  much  clarity  on  the 
figures  as  possible.       If,  in  future  applications  of  this  program,  these 
scales  are  not  convenient,   change  the  numbers  100  and  40  in  statements 
001620-001650  inclusive,  002910-002940  inclusive  and  004110-004180 
inclusive  to  more  appropriate  values. 

PART  9,  GRAPHICAL  OUTPUT 
Part  nine  of  the  program  consists  of  three  statements  which  call  the 
Subroutine  DRAW  for  plotting  the  experjunental  and  smoothed  values  of 
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Okj  around  the  circumference  of  each  pipe  section  at  a  particular  rosette,.i 
a  title  change  statement,  and  three  more  statements  which  call  the  Sub- 
routine DRAW  for  similar  plotting  of  T*y  values.  The  CALL  DRAW  statements 
occur  in  groups  of  three;  the  first  one  plots  a  circle  of  zero  stress, 
the  second  one  plots  the  smoothed  points,  and  the  third  one  plots  the 
experimental  values.  Since  compression  stress  appears  as  a  negative  stress 
while  tension  stress  appears  as  a  positive  stress,  a  polar  plot  could  not 
be  utilized  unless  the  position  of  zero  stress  was  far  enough  away  from 
the  center  of  the  plot  to  prevent  the  stresses  from  being  plotted  180 
degrees  out  of  position.  The  proper  scaling  of  the  plots  must  be  done 
before  the  program  is  utilized.  Appendix  D  ,  the  Subroutine  DRAW  writ eup, 
gives  full  details  on  both  scaling  and  what  should  be  placed  in  a  CALL 
DRAW  statement. 

PART  10,  COMPUTATION  OF  FORCE  AND  MOMENT  COMPONENTS 
Part  ten  of  the  program  computes  components  of  the  three  forces  and 
the  three  moments  at  each  ring  which  are  used  for  the  equilibrium  check. 
The  x-direction  is  taken  as  the  direction  of  flow  in  the  pipe.  The  y- 
direction  is  taken  as  the  direction  from  the  center  of  the  pipe  toward 
rosette  number  12.  The  z-direction  is  taken  as  the  direction  from  the 
center  of  the  pipe  toward  rosette  number  three.  XFOR  is  the  force  in  the 
x-direction  and  XMOM  is  the  moment  in  the  x-direction.  YFOR  is  the 
force  in  the  y-direction  and  IMOM  is  the  moment  in  the  y-direction. 
ZFOR  is  the  force  in  the  z-direction  and  ZMOM  is  the  moment  in  the  z- 
direction.  The  three  force  components  and  three  moment  components  are 
used  in  the  manual  computation  of  the  three  equilibrium  forces  and  three 
equilibrium  moments. 

PART  11,  GAUSS  2,  SOLUTION  OF  SIMULTANEOUS  EQUATIONS 
Part  eleven  of  the  program  is  the  Subroutine  GAUSS  2.  This  is  a 
standard  subroutine  in  the  computer  subroutine  library  at  the  U.  S.  Naval 
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Postgraduate  School.  The  computer  facility  explanation  of  this  sub- 
routine is  Appendix  C  of  this  thesis, 

PART- 12,  DRAW  *  GENERAL  PROGRAM  FOR  GRAPHICAL  OUTPUT 
Part  twelve  of  the  program  is  the  Subroutine  DRAW.  This  is  also  a 
standard  subroutine  in  the  computer  facility  subroutine  library  at  the 
U .  S.  Naval  Postgraduate  School.  The  computer  facility  explanation  of 
this  subroutine  is  Appendix  D  of  this  thesis. 
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4.  Data  Evaluation 

GAGE  ORIENTATION 

The  gage  orientation  is  one  of  the  important  items  to  which  careful 
thought  should  be  given  in  order  to  obtain  the  most  meaningful  data.  The 
most  important  consideration  concerning  gage  orientation  is  the  number  of 
rosettes  per  ring  and  the  spacing  between  rosettes.  As  previously  discussed 
in  section  2,  it  is  easier  to  obtain  Smoothed  values  of  CTn  and  t^f  when 
an  even  number  of  rosettes  equally  spaced  are  employed  because  equations 
with  only  one  unknown  are  more  readily  solved  than  several  equations  with 
several  unknowns.  The  actual  number  of  rosettes  per  ring  employed  in  the 
SSBN  640  instrumentation  appears  to  be  a  good  number  for  the  size  piping 
system  to  be  analyzed.  On  future  instrumentation  projects  the  size  of 
the  piping  system  should  be  considered  in  determining  the  even  number  of 
gage  rosettes  to  be  used. 

In  the  system  analyzed  the  rosette  orientation  was  shifted  45° 
alternately  between  even  and  odd  numbered  gages  (See  figure  1.2.1). 
This  shift  in  orientation  caused  minor  difficulties  in  obtaining  the 
proper  angle  between  the  maximum  principal  stress  and  direction  of  flow 
which  had  to  be  checked  very  closely;  it  is  also  a  source  of  possible 
misunderstanding.  The  identical  results  could  have  been  obtained  with  all 
rosettes  oriented  identically  and  thereby  eliminating  these  minor  diffi- 
culties. It  is  highly  recommended  for  future  instrumentations  that  all 
rosettes  be  oriented  identically. 

It  might  be  beneficial  to  consider  the  use  of  T-  DELTA  rosettes  for 
future  instrumentation  projects.  In  several  rosettes,  one  leg  indicated 
erroneous  readings.  If  a  fourth  leg  of  the  rosette  had  been  available, 
the  rosette  might  have  produced  meaningful  results  not  otherwise  obtain- 
able. 
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ACCURACY  CF  READINGS  -  DOCK  TRIAL  DATA 

The  readings  taken  during  dockside  hydrostatic  tests  were  recorded 

to  the  nearest  10yt^in/in.  Since  it  is  possible  to  read  an  SR-4  Indicator, 

Type  Nj  more  accurately  than  the  nearest  10  jL^Ln/in,  the  accuracy  of  all 

readings  is  not  considered  to  be  the  best  possible.  It  is  desirable  in 

future  tests  where  an  SR-4  Indi cat or^ Type  H} is  used  to  read  or  estimate 

strain  gage  indications  to  1  i(in/in.  Seven  pressure  increments  were 

recorded  during  the  tests. 

One  of  the  items  printed  out  in  Part  2  of  the  digital  computer 

program  is  DIFFS(I),  the  difference  between  each  actual  strain  gage 

reading  and  the  least  squares  calculated  best  line  fitting  valufcj.  In 

evaluating  the  actual  data,  any  value  of  DIFFS(I)  greater  than  10yi(Ln/in 

is  considered  to  indicate  an  inaccurate  reading.  With  this  as  a  basis 

of  inaccuracy,  the  following  statistical  tabulation  of  the  data  has  been 

made. 

Table  1 

Total  number  of  readings  taken  1939 

Total  number  of  readings  with  DIFFS(I)  greater  than  10 ^in/in  49 

%  of  total  in  error  2.52$ 

Total  number  of  gage  eleiaents  288 

Number  of  inoperative  elements  11 

Number  of  elements  with  one  reading  more  than  10/x;n/in  in  error  15 
Number  of  elements  with  two  readings  more  than  10  /Liin/in  in  error  9 
Number  of  elements  with  three  readings  more  than  10  /uin/in  in  error  4 
Number  of  elements  with  four  readings  more  than  10  JU±n/±ti  in  error  1 
Total  number  of  inaccurate  or  inoperative  elements  40 

%  of  total  elements  that  were  inaccurate  or  inoperative  13.89$ 

Total  number  of  rosettes  96 

Number  of  rosettes  with  at  least  one  inaccurate  reading  or 

an  inoperative  element  Ring  1  4 

Ring  2  3 

Ring  3  4 

Ring  4  3 

Ring  5  3 

Ring  6  7 

Ring  7  2 

Ring  8  k . 

Total  30 

%  of  total  rosettes  inaccurate  or  inoperative  31.25$ 

Total  number  of  rosettes  inoperative  due  to  faulty  element  9 
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$  of  total  rosettes  inoperative  due  to  faulty  element  9.375$ 

Ring  6  had  one  rosette  with  all  three  elements  faulty 

The  statistical  percentages  shown  in  Table  1  have  been  interpreted 

in  the  following  manners 

a.  Faulty  elements  caused  9o375$  of  the  rosettes  to  be  unusable 
in  evaluating  the  system*.  Although  thief  percentage  is  considered  to  be 
high  in  obtaining  good  experimental  results,  it  is  not  excessive  in 
evaluating  this  system. 

b.  The  percentage  of  inaccurate  rosettes  caused  by  a  reading 
error  of  10yUin/in,  21.875$*  does  not  have  any  real  meaning  because 
only  one,  two,  or  even  three  inaccurate  readings  might  not  have  caused 
the  slope  of  the  strain  depth  curve  to  be  enough  in  error  to  make  the 
rosette  give  a  faulty  overall  stress  calculation. 

c.  The  percentage  of  inaccurate  readings,  2.52$,  is  well  within 
the  experimental  accuracy  expected,  if  not  even  less  than  could  be  expected 
from  a  test  performed  under  the  conditions  of  this  test. 

PRINCIPAL  STRESS  RESULTS  -  DOCKSIDE  HYDROSTATIC  TESTS 
The  Bureau  of  Ships  Shipbuilding  Specif ication,  MILT  16420G,  states 
that  the  maximum  allowable  stress  in  70-30  copper  nickel  piping  under 
cyclic  loads  is  15000  psi.  In  publications  on  material  properties,  the 
^Leld  strength  of  70-30  copper  nickel  is  25,000  psi.  The  maximum  value 
of  (JJn^x  at  any  rosette  at  100$  test  depth  calculated  from  the  dockside 
hydrostatic  test  data  was  found  to  be  slightly  less  than  10,000  psi. 
This  indicates  that  the  maximum  principal  stress  does  not  reach  the  allow- 
able stress  and  is  well  below  the  yield  strength  at  any  instrumented  point. 
Since  the  rosette  ring  locations  were  chosen  to  be  representative  of  the 
entire  system,  this  is  sufficient  proof  that  in  all  probability  the 
stresses  within  the  piping  system  do  not  exceed  the  allowable  stress 
limit  of  the  material.  It  is  also  proof  that  the  stress  levels  were 
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sufficiently  low  to  assure  elastic  action,  as  is  assumed  in  the  analysis 
of  the  data. 

COMPARISON  OF  EXPERIMENTAL  DATA  TO  SMOOTHED  CURVES 

The  sixteen  figures  appearing  in  Appendix  E  are  polar  plots  of 
smoothed  axial  and  shearing  stress  around  the  circumference  of  the  piping 
at  each  rosette  ring  and  the  experimental  values  of  the  axial  and  shearing 
stress  at  the  individual  rosette  locations.  The  circle  ten  units  in 
radius  on  the  On  plots  and  four  units  in  radius  on  Txy  plots  is  the  circle 
of  zero  stress  on  these  plots.  For  the  CTn  plots  stresses  inside  the  sero 
stress  circle  are  Compression  while  those  outside  the  circles  are- 
tension*  .  For  the  UY  plots  stress  inside  the  zero  stress  circle  is 
shearing  stress  in  the  negative  direction  while  stress  outside  the  zero 

stress  circle  is  positive  shearing  stress.  The  scales  on  the  plots  are 
arbitrary  to  remain  within  the  unclassified  security  classification  as 
noted  in  section  1.  The  stresses  shown  are  plotted  per  percent  of  test. 
depth. 

It  is  readily  apparent  by  close  examination  of  the  sixteen  figures 
that  the  experimental  values  are  in  close  agreement  with  the  smoothed 
curves.  (The  experimental  points  are  marked  by  an  x  while  the  smoothed 
values  appear  as  a  closed  curve.)  The  x  at  the  center  of  some  plots  is 
to  be  ignored  as  it  is  only  a  necessary  step  in  having  the  CDC  1604 
computer  draw  the  plots.  It  will  also  be  noticed  that  some  plots  do  not 
have  all  twelve  experimental  points  shown.  The  missing  experimental 
points  are  an  indication  of  the  location  of  faulty  rosettes  (inoperative 
elements)  as  no  data  from  fully  operative  rosettes  was  rejected  in  this 
evaluation.  Had  the  missing  rosettes  been  operative,  the  shape  of  the 
curves  might  have  been  altered. 
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EQUILIBRIUM  RESULTS  FROM  DOCK-TRIAL  DATA 
Equilibrium  calculations  have  been  performed  on  the  segments  of  piping 

between  rosette  rings  1  and  2,  3  and  4,  and  7  and  8.  Calculations  were 
not  performed  on  segments  connecting  any  other  two  rings  because  unknown 
forces  due  to  foundations,  pipe  hangers,  or  other  constraints  exist 
between  these  segments.  Table  2  gives  a  tabulation  on  the  three  segments 
analyzed. 

Table  2 


(Resultant  moments  are  about  lower  numbered  end  of  segments) 


Theoretically,  for  equilibrium  to  be  satisfied,  all  values  of  Fr*., 
Ft^  Fe,»,  Mr*>  Miet,  and.  M^ should  be  zero.  Although  only  one  of  the 
eighteen  resultant  values  is  zero,  the  other  values  are  not  excessively 
high.  Any  one  of  the  following  influences .could  have  caused  the  non-zero 
values. 

a.  The  distances  between  rosette  rings  that  were  used  in 
calculating  the  resultant  moments  were  distances  measured  from  the  piping 
blueprints  using  the  approximate  locations  of  the  rosette  rings.  Since 
exact  values  of  these  distances  were  not  available,  these  were  the  best 
values  available  to  use  even  though  they  could  be  a  source  of  error  in  the 
calculations. 

b.  The  values  of  pipe  radius  and  thickness  used  in  calculating 
the  forces  and  moments  at  each  rosette  ring  were  average  values  calculated 
by  averaging  the  values  of  thickness  and  radius  at  each  rosette  within  the 
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ring.  In  addition  the  values  of  radius  at  each  rosette  reported  in 
Propulsion  Plant  Test  Fonn  721-5  were  obtained  by  measuring  pipe  diameters 
and  dividing  by  two,  which  is  not  a  true  value  of  the  radius  considering 
that  the  circumference  of  the  pipe  is  not  a  true  circle.  Since  outside 
diameter  measurements  for  rosette  ring  number  1  were  not  recorded,  the 
value  of  radius  of  ring  1  used  in  force  and  moment  calculations  was 
the  average  of  the  other  seven  rings. 

c.  The- value  of  the  forces  exerted  by  the  fluid  within  the  pipe 
on  the  inside  surface  of  the  pipe  was  determined  using  average  values  of 
radius  and  wall  thickness.  Using  average  values  instead  of  exact  values 
could  have  once  more  caused  errors. 

d.  Errors  in  the  Fourier  coefficients  caused  by  errors  in  strain 
gage  readings  could  also  have  contributed  to  errors  in  the  forces  and 
moments  calculated  at  each  rosette  ring. 

Considering  all  the  possible  sources  of  error  in  the  equilibrium 
calculations,  it  is  reasonable  to  expect  that  all  the  resultant  forces  and 
moments  would  not  be  zero.  However,  some  of  the  values  in  table  2  are 
much  larger  than  could  be  considered  reasonable  for  the  method  used. 
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5 •  Conclusion 

VALIDITY  CF  RESULTS 

As  indicated  by  the  results  shown  in  this  thesis,  the  method 
developed  can  readily  be  used  to  obtain  stress  information  concerning 
the  system.  However,  the  validity  of  the  equilibrium  results  will  depend 
largely  on  the  accuracy  of  the  strain  gage  data  and  other  data  needed  for 
the  equilibrium  calculations. 

In  order  to  increase  the  accuracy  of  the  data  and  thereby  increase 
the  validity  of  the  results,  the  following  recommendations  are  submitted: 

a.  Appoint  a  project  engineer  with  experience  in  piping  theory 
and  strain  gage  measurements  to  have  overall  control  of  the  instrumentation 
of  the  system,  data  taking  during  all  tests,  and  reporting  of  results 
obtained  from  the  data.  This  one  person  would  supervise  the  project  from 
the  start  of  instrumentation  to  the  submission  of  the  final  report.  Having 
one  man  in  charge  of  the  project  from  start  to  finish  will  greatly  reduce 
the  chances  of  faulty  communication  between  people  working  on  the  project, 
reduce  the  chances  of  necessary  data  items  being  omitted  from  the  data 
taken,  and  reduce  the  chances  for  error  in  the  actual  data  taken. 

b.  Utilize  automatic  scanning  and  recording  equipment  to  record  all 
strain  gage  data  during  actual  tests.  Automatic  recording  equipment  can 
observe  strain  gage  reading  much  more  accurately  than  equipment  read 
manually  and  in  much  less  time.  In  addition,  human  recording  errors  could 
be  eliminated.  A  rapid  record  of  strain  gage  reading  at  any  depth  could 
also  be  taken  using  the  automatic  scanning  and  recording  equipment. 

c.  An  important  aid  to  the  person  analyzing  the  experimental 
data  from  these  tests  would  be  notes  taken  during  the  period  that  data 
is  being  taken.  Some  of  the  information  that  would  be  useful  as  notes 
is  any  unusual  conditions  that  occurred  during  data  taking,  comments 
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concerning  any  data  that  might  give  an  immediate  indication  as  to  being 
inaccurate  or  unreliable,  the  conditions  that  existed  during  the  periods 
between  the  taking  of  data,,  etc, »  The  project  engineer  should  be 
observing  all  tests  to  take  his  own  notes  as  well  as  having  notes  taken 
by  all  technicians  involved  in  the  tests, 

d.  Additional  data  to  be  taken  during  the  tests  should  include 
actual  measured  distances  between  two  adjacent  rosette  rings  for  use  in 
more  accurate  equilibrium  calculations  and  the  location  of  the  welled 
seams  for  each  section  of  pipe. 

By  further  refinement  of  data  taking  methods  and  actual  data  taken, 
the  method  of  evaluation  developed  in  this  thesis  can  be  used  as  a  tool 
in  analyzing  any  piping  system  instrumented  in  a  manner  similiar  to  the 
Main  Seawater  Piping  System  on  SSBN  640. 

SUGGESTED  FUTURE  REFINEMENTS  FOR  METHOD  OF  EVALUATION. 

-  •  ....  mm- 

Certain  refinements  to  the  present  method  of  evaluation  have  been 
considered;  however,  sufficient  time  is  not  .available  to  permit  these 
ideas  to  be  included  in  this  thesis.  The  basic  content  of  these  ideas 
is  listed  in  the  following  paragraphs. 

Part  2  of  the  program  could  be  expanded  to  include  a  method  which 
would  investigate  the  value  of  DIFFS(I),  This  would  be  a  two  stage 
process.  The  zero  percent  depth  intercept,  aj  slope,  bj  and  values  of 
DIFFS(I)  would  be  found  as  they  are  in  the  present  program.  If  DIFFS(I) 
were  greater  than  a  predetermined  limit  value,  the  data  which  produced 
that  DIFFS(I)  value  would  be  rejected  and  new  values  of  aero  percent  depth 
intercept  and  slope  would  be  found.  An  addition  of  the  type  proposed 
would  improve  the  accuracy  of  all  resultant  calculations. 

The  present  theory  might  be  further  developed  to  find  a  method  to 
make  some  use  of  strain  gage  rosettes  which  have  one  or  two  inoperative 
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elements.  When  the  center  element  is  the  only  one  inoperative,  some  useful 
information  may  obviously  be  obtained  from  the  rosette.  If  the  orientation 
is  that  of  the  even  numbered  rosettes  as  shown  in  figure  1.2.1,  elements 
A  and  C  operative  and  B  inoperative  would  permit  evaluating  the  axial 
stress.  If  the  orientation  is  that  of  the  odd  numbered  rosettes  as  shown 
in  figure  1.2.1,  elements  B  and  C  operative  and  A  inoperative  would  permit 
evaluating  a  shear  strfess.  Beyond  this,  if  rosette  failure  were  other 
than  as  indicated  above)  there  is  probably  some  statistical  way  of  making 
use  of  the  data  afforded  by  the  operative  elements.  However,  the  theory 
of  such  usage  remains  to  be  worked  out. 

These  two  refinements  added  to  the  analysis  method  might  produce  more 
meaningful  results  from  the  dockside  hydrostatic  test  data. 
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APPiiNDIX  A 
COMPLETE  PROGRAM  LISTING 
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,.  JO^ri&S-ar,   <CCH»  J  W 

program  strain 

^ I  MANSION  X  (  8  )  ,  X^-vyF  !  3  )  »Y  (  8  )  »  YSA  vr  (  ?.  I  ♦  v-  1  T  (  3  j  ,')\  czz  [ 
1  .  r-.-'  ( 1 1 ) »CK ( 1 1 ) » A ( 1 1 ♦ 1 1 >  »3N ( 1  1„  1 ) ,CN ( 1 1 » 1 ) »ZX ( 1 1 , 1  > , X 
2),~HI1{63C)»SIGXX(12)  ,TA'JXY(  i  ;  )  ,  I  T  !  t  Lr  ,  :  ,  ;  ,  S  I  GN7  !  £  3  " 
3MAX(  12) »SMIN< 12) »4NG! 12) »THETA( 12) »XSXX( 12)»VSXX( 12) 
4  (  12  )  »XSNR  (  :>3'N  !  ,  YS\R(  £3:  )  »XSN7  '  630  )  »VSN?  (  6  30  )  »X  r.\R<  53 


.  ?(  3) 
11.1) 


*  v  (  12 


»TH(S) ■ 
P  A  R  T  1 , 


.(8)  , 

PAR' 


XFORf 3 ) » YFORI 8 

i  rea; 


T  ' 


5T."  M630)  »YT,\7(63V') 
6V(8  )  »Z  ''GM(  9  ) 
ThlT  !  c-  THE  ST-RT  C" 
.-  r  A  D  3  9  5  ♦  N  »  E  » U 
999  FOR,'-' A;  (  114, 2F  15.5) 
TH( 1 1=0.625 
."••  A  (  1  )  =  6  .  4  r 
TH( 2  1=0.634 
RA(2)=6.390 
THt  3  )=0.620 
RAi  3 1=6.409 
THt  4) =0.633 
RA(4)=6.374 
TH( 5 1=0.428 
RA(5 1=6.416 
THt  b  1  =0.64  0 
RA (6  1=6.354 
TH( 71=0.423 
RA{ 7)=6.411 
TH(8)=0.423 
RA(8  1=6.410 
DO  6  5  K  =  1  »  8 
DO  3  1=1,28 

3  AM(  I  1  =  0. 
DO  4  1=1,7 
BM( 11=0. 

4  CM( 1 )=0. 

Z  =  3. 14159/6. 
THIS  IS  THE  >^D  OF  PART  1. 
13  DO  15  L=l,12 

THIS  IS  THE  START  OF  PART  2.   PART  2  READS  THE 
PROGRAM  TO  DETERMINE  THE  STRAIN  PFR  PERCENT  OF 
INDIVIDUAL  GAGE. 
57  DO  58  M=  1,3 

READ  100,  (X(I),Y(I),I=1,N) 
10.   FORMAT ( 14F5.0) 

DO  996  I=1,N 

XSAVE1 I )=X( I ) 
99b  YSAVF (  I  )=Y<  I ) 

YSUM    =    0. 

SUMX    =    0. 

SIJMY     =     0. 

r,i[.'n  =  D. 
SlJMXX  =  (J. 
XN    =    0 

DO    30    I    =    1,    N 
SljMX  =  SUMX+X  (  I  1 
,     ■■  1    -    SUMY+Y  (  I  ) 

;xx  =  suMXX  +  jx(  i  )*Xi  i » j 


,  c .  ;_ 


\  ■•■  I 


-t  ) 
(  12 
1  ,S 

T  ^   V 

)  ,x 


DATA  AT  THE  END  OF  THF 
DEPTH  FOR  EACH 


"",  /\  *>  f*  *    *> 

1 0 1  n  2  0 

3  0  0  0  3  C! 
n  o  n  c  ^  o 

"C0C6C 

CCO 070 

ooooso 

0  0  0  0  9  0 
COO  100 
000  HO 
090  120 
POOl 30 
000  1 40 
000  150 
000160 
000170 
000180 
000  19C 
000200 
000210 
000220 
000230 
000240 
00^250 
000260 
000270 
000280 
000290 
000300 
000310 
000320 
000330 
000340 
000350 
000360 
000370 
000380 
000390 
000400 
000410 
000420 
000430 
000440 
000450 
000460 
000470 
000480 
000490 
000500 
000510 
000520 
000530 
000540 
C00550 


46 


30  SUMXY=SUMXY+(X( I )*Y( I  )  )  000560 

DENOM  =  XN*SUMXX-SUMX*SUMX  000570 

AHUM  =  SUMY  *  SUMXX  -  S'JMX  *  SIJVXY  000580 

AZ  =  ANUM  /  DFNOM  000590 

RNl.iM=  XN*SUMXY-SUMX*SUMY  000600 

R(M)  =  BNUM  /  DENOM  000610 

DO  40  J  =  l.N  000620 

YC  =  AZ  +  B(M)  »X(J)  000630 

DIFF  =  Y(J)  -  YC  000640 

40  YSUM  =  YSUM  +  DIFF*DIFF  000650 

S2  =  YSUM  /  (XN  -  1.)  000660 

SIG=S2/DENOM  000670 

SRA=SQRTF ( SUMXX*S I G )  000680 

SRB=SQRTF(XN*SIG)  000690 

PRINT  997,  AZ,  R(M),  SRA,  SRB  000700 

997   FORMAT (24HC   STRAIGHT  LINE  FITTING//  25H  USING  THE  FORM  Y=BX  +  A,    000710 

1/  4H  A  =  F14.8/  7H  RIM)  =  F14.fi/  26H  STANDARD  DFVIATION  OF  A  =  E14   000720 

2.8/  29H  STANDARD  DEVIATION  OF  R(M)  =  E14.8)  000730 

801  SUMFIT=.0  000740 

DO  90  1=1, N  000750 

YC1T(I)  =  B(M)*X(I)+AZ  000760 

D1FFSI  I  )=YFIT  I  I  )-Y( I  )  000770 

90  SUMFIT=SUMFIT+  DIFFS(I)  000780 

995  PRINT  99,SUMFIT,(X( I  ) ,YSAVE( I  )  ,YI I  )  .YFITI  I  ),DIFFS( I ) ,1  =  1, N)  000790 

99  FORMAT (9H0SUMF IT  = E 1 2 . 5// / 10X1HX 14X 1HY9X 10HYTRANSPOSE8X4HY F I Tl 1 X5H   000800 

1DIFFS///I 5E15.5) )  000810 

58  CONTINUE  000820 

C   THIS  IS  THE  END  OF  PART  2  AND  THF  START  OF  PART  3.   THIS  PART  000830 

C   DETERMINES  THE  PRINCIPAL  STRESSES  AND  THE  ANGLE  FROM  THE  PIPE  000840 

C   CENTERLINE  OF  THE  MAXIMUM  PRINCIPAL  STRESS  FOR  EACH  ROSETTE.  000850 

IF(BID)  401,410,401  000860 

401  IFI3I2))  402,410,402  000870 

402  IFIBI3))  403,410,403  00C880 

403  AA=E/2.  000890 
BA=B(1)  +  8(3)  000900 
CA=B(1)  -  8(3)  000910 
D=l.+U  000920 
F=l.-U  000930 
G=2.*B(2)  -  BA  000940 
H=DA/F  000950 
P=1./D*SQRTF( CA*CA+G*G)  000960 
SMAX(L)=AA*(H+P)  000970 
SM!N(L)=AA*(H-P)  000980 
ANG(L) =0.5*1 ATANF(G/CA)+( 3.14159/2. )*( 1. -I CA/A3SF(CA))))*180./3. 14   0009  90 

1159  001000 

GO  TO  15  001010 

410  SMAX(L)=0.  001020 

SMIN(L)=0.  001030 

ANG(L)=C.  001040 

15  CONTINUE  001050 

DO  115  L=l , 11 ,2  001060 

115  THETA(L)=ANG(L)-45.  001070 
DO  116  L=2,12,2  001080 

116  THF.TA(L)  =  AN3(L)  001090 
PRINT  5, ( SMAX (L ) , SMINI L ) ,THLTA ( L) ,L=i ,12 )  001100 

5  FORMAT! ///6X7HSMAX (L) 8X7HSMIN( L ) 8X8HTHETA ( L) /// ( 3 El 5. 5 )  )  001  110 
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AND  Trie  START  OF  ^ART 
SOLUTION  TO  SOLVE  T HE 


)  )*COSF(THETA(L  )  ) 
*0.5*SINF(2.*THETA(L 


THI  S  IS  THE  END  OF  r-ART  3 
WHETHER  THE  SHORT  OR  LO'vG 
MUST  3E  USED, 

DO  420  L=ltl2 

IF (SMAX(L) )  42C»t3?t <+20 
42r  CONTINUE 
THIS  IS  THF  END  OF  PA^T  4  AND  THF  START  OF  PART  5 
TH^  VALUES  OF  THE  COEFFICIENTS  US^D  IN  THE  THEORE 
OF  NORMAL  STRESS  AND  SHEAR  STRESS  BY  THE  SHORT  OR 
GAGES  EQUALLY  SPACED  METHOD. 

DO  16  L=1.12 

THETA(L)=THETA{L)*3.14159/18C- 

SIGXXt  L)  =  SWAX(L)«COSF( THETA1L 
ID  )*SINF(THETA(LI  ) 

TAJXYI L)=( SKAXI L ) - SM I N ( L 

XL  =  L 

PHI (L)=XL*Z 

SI=SINF(PHI  (L)  ) 

CI=COSF(PHI  (L)  ) 

SI2  =  SINF( 2.*PHI (L)  ) 

CI2=COSK(2.*PHI  (L)  ) 

SI3=SINF( 3.*PHI (L) ) 

CI3=COSF( 3.*PHI (L)  ) 

SISQ=SI*SI 

CISQ=CI*CI 

S2SQ=SI2*SI2 

C2SQ=CI2*CI2 

S3SQ=SI3*SI3 

C3SQ=CI3*CI3 

AM(1 )=12. 

AM(6)=AM( 6)+SISQ 

AM(  10)=AM(  1O+CIS0 

AM( 13)=AM( 13)+S2SO 

AM(  15)=AM( 15)+C2SO 

AM(26)=AM(26)+S3SQ 

AM(28)=AM(28)+C3SQ 

BM(  1  )  =  BM(  D+SIGXXI  L) 

BM(2)=BM( 2)+SI*SIGXX(L  ) 

BM(3)=BM(3)+CI*SIGXX(L ) 

BM(4)=BM(4)+SI2*SIGXX(L) 

BM(5)=BM(5)+CI2*SIGXX(L) 

BM(6)=BM(6)+SI3*SIGXX(L) 

BM(7)=BM( 7)+CI3*SIGXX(L> 

CM( 1 )=CM(  1)+TAUXY(L) 

CM(2)=CM(2)+SI*TAUXY(L) 

CM(3)=CM(3)+CI*TAUXY(L) 

CM(4)=CM(4)+SI2*TAUXY(L) 

CM(5)=CM( 5)+CI2*TAUXY(L) 

CM(6)=CM(6)+SI3*TAUXY(L ) 

CM(7)=CM( 7)+CI3*TAUXY(L> 

XSXX(L)=( 100.  +  SIGXX(L>  >*CI 

YSXX(L)=( 100.+SIGXXI L> )*SI 
XTXY(L) =(40.+TAUXY(L) ) *C I 

YTXY<L)=(40.+TAUXY(L) I *S I 
16  CONTINUF 

PRINT  7.(SIGXX(L).TAUXY(L) »L=1»12) 


PAf 
TRESis 


T  4  DETE 
DISTRIdU 


.  °ART 

TICAL 

EVEN 


5  DETER 
DETERMIN 

NUMBER  0 


RMINE; 

TION 


VINES 
AT  I  ON 
F 


+SMINI L )*SINF{ ThETA( 
)  ) 


001120 
001  130 
001140 
001150 
001160 
00117C 
001  130 
001 190 
001200 
001210 
C01220 
001230 
001240 
001250 
001260 
001270 
001280 
001290 
001300 
001310 
001320 
001330 
001340 
001350 
001360 
001370 
001380 
001390 
001400 
001410 
001420 
001430 
001440 
001450 
001460 
001470 
001480 
001490 
001500 
001510 
001520 
001530 
001540 
001550 
001560 
001570 
001580 
001590 
001600 
001610 
001620 
001630 
001640 
001650 
001660 
001670 


4d 


FORMAT ( / / /6v3hc 

•'  ( 1 . i 1  =  a  '•■  ( ;  i 

M2»2l=AM(6) 

A  (  3  *  3 ) =  AV ( 1 C ) 
A ( 4  .  4 ) =  AX { 1 3 ) 

A  (  5  ,  5  )  =  A  ■  '■  (15) 
A  (  6  » 6  )  =  A  V  (  2  6  ) 

a:  "',?)  =a:-m  28) 


I  r.  v  y  i  ">  y  i-  h  ""  A  i  v  v  /  /  /  f 


i  ) 


: ) 


DO  11 

r  n  ( ; » ] 
cm 1 1  ,  i 

2  X  (  1  ,  ] 


2,1 


2X( 

2 X  I  3  ,  i 
2X14,1 
2  X  (  5  «  1 
2X (6,^ 
2X( 7,i 
Xv (  1  ,1 
Xy (2.1 
XY ( 3,  1 
X  Y  (  4  ,  1 
X  Y  (  5  »  i 
XY (6,  1 
XY<7,1 

pr:nt 


■  i 


.96 


156 


TH 
TH 
OF 
DE 
UN 
430 

42 


4h  1 


=  ?y<  :  i 
=cm<  ; ) 
=  BN(  l  ♦: ) 

=  o  \  (  2  »  1  ) 

=  6  N  •!  3  » 1  ) 

-  b  N  i  h  ,  1  ) 
=  5N(  5.1  ) 
=bN( 6.1 ) 
=  t3  N  (  I' .  i  ) 

=cn:  i.i i 

=  CN( 2.1 ) 
=  C  N  (  3  »  1 ) 
=  C  N  (  <t  » 1  ) 
=CN( 5»1) 
=CN( 6.1 ) 
=CN( 7»1) 
196.2X1 1 , 

FORMAT I///11H  2 
iE14.8//llH  Z X  (  4 
2.8//11H  2X(7.i) 

PRINT  li6*XY(li 

FORMAT ( ///11H  X 
lE14i8//ll'H  XYC4 
2.8//UH  XYI7.1) 

GO  TO  406 
rS  [S  THF  FN~>  OF 
r:  VALUES  OF  THE 

THE  VALUES  OF  T 
ThRMlNATION  OF  N 
EQUALLY  SPACED  G 

00  46  L=i»12 

IF (SMAXfL ) )  44^ 

S IGXX(L)- 

TAUXY(L)= 

XSXX(L) =0 

YSX  X ( L) =  0 

XTXY(L)=0 

YTXY(L)  =0. 

GO  TO  4  6 

THFTA( L)=THFTA( 

SIGXXt L >=SMAX< L 
ID  )*SINF  (THFTA( 

TAUXY (L)= (SMAX( 

XL  =  L 

PHI (L)=XL>2 

SI=S1NF (PHI (L) ) 

CI=C05F(PHI  (L  )  ) 


.4) 


.  - 

♦  3 

.4 

,5 
.  6 

♦  < 
X  ! 
1  ) 
= 

14 
Y  ( 
i  ) 
—     i 
14 


l)*ZX(3fl)(Z 
=  E14.8//11H 

14.6//11H  Z X ( 
8//  ) 

.  1  i  »  X  v  i  3  ,  i  )  ,  x 
-  E14.3//11H 
14.8//11H  XY( 
8  /  / ) 


PAR'T  5  AND  THF  STAR 
COEFFICIENTS  USED  IN 
HE  COEFFICIENTS  USED 
OhMAL  STRf SS  AND  SHE 
AGES  M£THOl). 

,42  ,440 


X  <  4, 
ZX(2 
5.  i  ) 

Y(4» 

XY(2 
5.1) 


T  OF 
7HE 
IN 

Ai\  S 


i:  .2 
,1 ) 


l )  »x 

« i ) 


PAR 
THE 

THE 

TRES 


x  ( ■;> .  i 
=  Ei4 

14.6/ 

Y(  5»1 

=  El  A 

14.3/ 


6. 
OKETI 

TrlFO!-' 
S  BY 


)  .  Z  X  f  6  . 

. 0/ / lln 

/llri  ZX 

/  *  X  Y  (  6  t 

.0//11H 
/11H  XY 


PART  6 
CAL  DET 
ETiCAL 
THE  LON 


1).Z 

2X( 

(6.1 

1)  .X 

XY( 

(6.1 


X.  (  7  ,  1  ) 

3,1  )  = 
)  -  El  4 

Y (7,  1  ) 
3,1  )  = 
)  -     E14 


DET 
ERMI 

G  OR 


ERMINE' 

NAT  ION 


D*  3.  14159/180, 

)*COSF ( THETA ( L ) ) *COS 

L)  ) 

D-S  'IIN(L)  )  *0.5*C'INF 


F  (  T ,  I 
(  2  .  •' 


E  T  A  < 

InLT 


L)  )  *■' 
A(l) 


iIN(L').*SINFt  ThETAI 


0  0 1 6  8  n 

0  3  1  6  -  0 
001700 

301 710 

C  C  i  7  2  0 
331 7  30 
OOi740 
301 7  50 
SCI  763 
: S17  70 
0  017  8". 
-)01 7  on 
001900 
001810 
301823 
001630 
301840 
001850 
001860 
001870 
001883 
^01390 
CO  1900 
0  019  10 
001920 
0010  30 
00i940 
0019  50 
001960 
001970 
001980 
0  019  90 
002000 
002010 
002  020 
002030 
002040 
002050 
002  060 
302070 
002080 
002C90 
002100 
002110 
0C2120 
00  2  130 
002  140 
002  150 
002160 
002170 
002180 
00219C 
00<:20G 
002210 
002220 
002230 
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". !  2  =  s ; 

3:3=3; 

C  t  3  =  r  c 
3! SQ  =  3 
CISG=C 
SC=SI* 
SIS2=S 
SIC2=S 
CIS2=C 
CIC2=C 
S2SQ=S 
S2C2=S 
C2SQ=C 
S I  S  3  =  S 
S I C  3  =  S 

crs3=c 

;.  i  c  3  =  c 

52S3=S 

52C3=S 
C2S3=C 
C2C3=C 
S3SQ=S 
S3C3=S 
C3S0=G 
AMI  1  ) 
AM (2  )  = 
AM(3)= 
AM ( 4 ) = 
AM( 5)= 
AMI  6)  = 
AMI  7)  = 
AM (8  )  = 
AM  19)  = 
AMI  10 
AM!  11 
AMI  12 
AM  (  i  3 
AM(  14 
AM  I  1  5 
AMI  16 
AMI  17 
AMI  16 
AMI  19 
AMI20 
AMI21 
AMI  22 
AMI  23 
AM  (  24 
AMI  2  5 
AMI  26 
AM  (  2  7 
AMI28 
BMI 1 )= 
BMI  2 )= 
BMI3)= 


F-12. 
F  I  2  . 

n  t  :  3 . 

5  F  (  3  . 
I*5i 
I*C1 
C! 

I  ■>■  £  I  2 
i  *  C  1 2 
1*5  12 
I*C12 
12*51 
I  2*  CI 
I2*CI 
I  *  S  i  3 
!  *  C  !  3 
I  *  S  I  3 
I --CI  3 
I  2*SI 
I  2  *  C I 
I  2*Si 
I  2 -CI 
I3*SI 
I3*CI 
I3*CI 
AMI  1  ) 
AMI  2) 
AM  I  3  ) 
AM(4) 
AMI  5) 
AMI  6) 
AMI  7) 
AM(  3) 
AMI  9) 
=  AM(  i 
=  AM(  1 
=  AM(  1 
=  AM<  1 
=  AM|  1 
=  AM(  1 
=  AM( 
=AMI  1 
=  AMI  i 
=AM(  1 
=  AM(  2 
=  AM(  2 
=  AM(  2 
=AMI2 
=  AM(2 
=  AM(2 

=  AM(  2 
=  AM(2 
=  AM<2 
BMI  1) 
BMI  2) 
BMI  3) 


P 

Oh  l 


!L)  ) 
I  L  )  I 
(Li  ) 


i 
3 
3 

'3 

3 

3 

3 

+  1. 

+  SI 

+  CI 

+  SI2 

+  CI2 

+  SI5Q 

+  5C 

+  SIS2 

+SIC2 

0)+CISQ 

11+CIS2 

2 )+CIC2 

31+52SO 

4)+S2C2 

5)+C2SQ 

61+SI3 

71+CI3 

3)+SIS3 

9)+SIC3 

01+CIS3 

1 )+CIC3 

21+S2S3 

3)+S2C3 

4J+C2S3 

5I+C2C3 

6)+S3rQ 

7  1+S3C3 

a >+C3SQ 

+  .S1GXXI  L) 

+  51 *SIGXX< L  ) 

+CI*5IGXX( L  ) 


0C22«+0 

*^2  2^ 

002260 

0  0  2  2  ■  "■• 

3  32260 

002290 

002 3CC 

30^310 

032320 

0C2330 

002340 

002350 

002  360 

002  370 

002380 

002  390 

002400 

002410 

002420 

002430 

002440 

002450 

002460 

002470 

CC2480 

002490 

002500 

002510 

002520 

002530 

002540 

002550 

002560 

C02570 

002580 

002590 

002600 

002610 

002620 

002630 

002640 

002650 

002660 

002670 

002680 

002690 

002700 

002710 

002720 

002730 

002740 

002750 

002760 

002770 

002780 

002790 

50 
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7  / 


""•  |  5 

-,:•(  7 

'■'(  I 
."•'(,.' 

-"(  ? 

C  M  (  U 
CM(  5 
CM  (6 
CM  I  7 
XSXX 

Y<":XX 

XTXY 
YTXY 
CONT 
PPIN 
FORM 

A(  1 

A(  1 
A(  i 

A(  X 

A(  1 
A(  1 
A(  1 
A(  2 
A<2 
A(2 
4(2 
A<2 
A(2 
A(2 
A(  3 
A(  3 
A(  3 
A(  3 
A(  3 
At  3 
At  3 
A(4 
A!4 
A14 
A(4 
A(4 
A(4 
A(4 
At  5 
Mb 
At  5 
At  5 
At  5 
At  5 
At  5 
A(6 
A(6 
A(6 


'.  (  u  )  +  5 1  2  * 
(51+CI2* 

(  A  )  +  ^  ;  ",  i 

'(  7)  +CI^* 


IGX  X ( L ) 

igxxi  l  ) 

I 'XX ( L  ) 
,  IC-XXt  L  ) 


rw( 1)+TAUXY(L) 

Cf.'t  2  )  +3  I* 7  A  :XY(Li 

r\'( 3 ) +c ;»:; ;yv ( L ) 

CM(4)+S I 2<Ti  JXY(L ) 

CM(  5)+CI2*TVJXv ! L ) 

CM( 6) +31 3*TA JXY(L) 

CN'.m+CI  3*TAUXY(L  ) 

)  =  (  1Q0.+S1GXX4  L)  )*CI 

)=(10O,+SIGXX{Ll >*SI 

)  =(  40.+TA'.:xv  (L)  )*CI 

)  =(  40  .+TA'.IXY  (  L)  )  *c  I 

'JF 

77 ,  (  5 1 GXX  (  L )  »T  AIJXY  (  L  )  »  L=  1  ♦  1 2  ) 

(  ///6X5nc-IGXXl'.'X5HTA'JXY///(  2E15.5 )  ) 

=  A  M  (  1  ) 

=AM(2 ) 

=  AM( 3  ) 

sAM(4) 

=AM(5J 

=AM< i6) 

=  AV.  (  17) 

=  Av,(2  ) 

=AM(fe) 

=AM(7) 

=AM{ 3) 

=AV(9) 

=  AM( 18) 

=AM( 19) 

=AM( 3) 

=AM(7) 

=AM( 10) 

=AM( 11) 

=  AM( 12) 

=AM(20) 

=  AV(2l  ) 

=  AM(4  ) 

=AM(8l 

=  AM!  11) 

=  AMI  13) 

=AM( 14) 

=AM(22) 

=AM(23  ) 

*AM(5 ) 

=AM(9) 

=AM( 12) 

=AM( 14) 

=  AM(  15) 

=AM(24 ) 

»  A  M  (  2  5  ) 

=AM( 16) 

=AM( 18) 

=AM(20) 


002600 
302810 
■)?2~2" 
^02830 
0O28-O 
002850 
002360 
00287n 
0C2880 
002890 
002900 
0029 i0 
00  2920 
002930 
002^40 
002950 
002  960 
002970 
002980 
002990 
003000 
003010 
003020 
003030 
003040 
003050 
00  3  060 
003070 
003080 
003090 
003100 
003110 
003120 
0C3130 
003140 
003150 
003160 
003170 
003180 
00  3190 
003200 
003210 
003220 
003230 
003240 
003250 
003260 
003270 
003280 
003290 
003300 
003310 
003320 
003330 
003340 
003350 
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A(6»4)  =  AM  (  2  2  ) 

A (6,5) a AMI  24) 

A(6.6)=AM(26) 

A  (  6  ,  7  )  =  A  X  (  2  7  ) 

A  (  7  •  I )  =  A  M  (  i  7  ) 

A( /,2) =AM< 19) 

A(  7.3)=AM(21  ) 

A  (  7  ,  4 )  a  A  M  {  2  3  ) 

A(7»5)=AM(25) 

A( 7,6)= AMI  27) 

A(7»7)=AM(28) 

PRINT  299 »( (A(J,I ) ,1=1 »7) »J 
299  FORMAT ( ///2X6HAI 1,1 )4X6HAi I 
1  I  »  6  )  4X6HA (I,7)///I7E10.4)) 

DO  i 17  1=1,7 

BN( 1,1 ) =  BM(  I  ) 
117  OX  1  ,1  )=CM(  I  ) 

NN  =  7 

MM  =  1 

EP=l.F-9 

CALL  GAUSS2 I NN , MM , EP » A ,BN » Z 

PRINT  799,ZX(1,1) »ZX(2»1) »Z 
799  FORMAT) ///11H  ZX(i.l)  =  E14 
1E14.8//11H  ZX(4,1)  =  E14.8/ 
2.8//11H  ZX(7,1)  =  E14.8//1 

CALL  GAUSS2 I  NN,MM, EP »A ,CN ,X 

PRINT  759»XY( 1,1) ,XY(2»1) ,X 
7  59  FORMAT ( ///11H  XY ( 1  ,  1)  =  El 4 
1E14.8//11H  XY(4,1)  =  E14.8/ 
2.8//11H  XY(7,1)  =  E14.8//) 
THIS  IS  THE  END  OF  PART  6  AND 
THF  TITLES  THAT  ARE  PUT 
456  DO  107  1=1,12 
107  ITITLE(I)=8H 

ITITLEI 1 )=8H  KOCH,  J 

ITITLEI 2) =8H.  W.»  JR 

ITITLE(7)=8H  SIG  VS. 

ITITLE(8)=8H  PHI 

ITITLEI 9)=8HRING  NO. 

GO  TO  (  301,302,3i~!3,3n4,305, 

301  ITITLEI  10)=8H  ONE 
GO  TO  315 

302  ITITLE(10)=8H  TWO 
GO  TO  315 

303  ITITLEI 10)=8H  THREE 
GO  TO  315 

304  I  TITLE!  10)=8H  FOUR 
GO  TO  315 

305  ITITLEI 10)=8H  FIVE 
GO  TO  315 

306  ITITLEI 10)=8H  SIX 
GO  TO  315 

307  ITITLEI 10)=8H  SEVEN 
GO  TO  315 

308  ITITLEI !0)=8H  EIGHT 
3] 5  ITITLFI  11  )  =  8H  3/7/66 


=  1,7) 

»2 ) 4*bHA 


I  ,  3  )4X6riA(  i»^;4X6riA(  I»5)4X6mA( 


X,K^R) 

X(3,1)»ZX(4,1),ZX(5,1),Z 
.3//11H  ZX(2,1)  =  Fi4.8/ 
/ 1 1 H  Z X ( 5 , 1 )  =  E14.8//11 

Y.KER) 

Y<3»1)»XY(4»1)»XY<5»1)»X 
.8//11H  XY(2,1)  =  E14.8/ 
/11H  XY(5,1)  =  E14.8//11 

THF  START  OF  PART  7.   ?A 
ON  THE  PLOTS  DRAWN  3Y  THF.  PROG 


X  (  6  ,  1 )  »  Z 

/11H  ZX ( 
H  ZX(6,i 


Y(6,l) ,X 
/11H  XY  I 
H  XYI6, i 


X ( 7,  1  ) 
3,1  )  = 

)  =  E14 


Y(7,l) 
3,1)  = 
)  =  E14 


RT  7  READS  IN 
RAM. 


306,307,3^8) ,< 


00336C 
0C3370 
003360 
003390 
003400 
00341C 
003420 
003430 
003*40 
0C3450 
003  460 
003470 
003480 
003490 
003500 
003510 
003520 
00  3  5  30 
00  3  540 
003550 
00  3  560 
003570 
003580 
003590 
003600 
003610 
003620 
003630 
003640 
003650 
003660 
003670 
003680 
003690 
003700 
003710 
003720 
003730 
003740 
003750 
003760 
003770 
003780 
003790 
003800 
003810 
003820 
003830 
003840 
003850 
003860 
003870 
003880 
003890 
003900 
003910 
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LB  =  4H  003920 

LD  =  8H  003930 

C   THIS  IS  THE  END  OF  PART  7  AND  THE  START  OF  PART  8.   PART  8  COMPUTES  003940 

C   THE  THEORETICAL  VALUES  OF  NORMAL  STRESS  AND  SHEAR  STRESS  AROUND  THE  003950 

C   PIPE  DIAMETER  AND  TRANSFORFS  THEM  INTO  VALUES  THAT  ARE  TO  BE  PLOTTED.  003960 

ZZZ=6. 28318/628.  003970 

DO  26  1=1,628  003980 

Z I = I  003990 

PHIKI  )=ZI*ZZZ  004000 

S0  =  SINF(PHI1( I )  )  004010 

C0  =  C0SF(PHI1<  I  ) )  004020 

S02=SINF(2.*PHI1( I  )  )  004030 

C02=C0SF(2.*PHI1( I  )  )  004040 

S03=SINF(3.*PHI1( I ))  004050 

C03=COSF(3.*PHI 1( I ) )  004060 

SIGN7( I )=ZX(1»1)+ZX( 2»l)*SO+ZX(3»l)*CO+ZX(4,l)*S02+ZX(5,l)*C02+ZX(  004  070 

16,1)*S03+ZX(7,1)*C03  004080 

TAUN7I I )=XY( 1,1 )+XY( 2»1 )*SO+XY(3»l )*C0+XY(4»1 )*S02+XY(5»1 )*C02+XY<  004090 

16,1 )*S03+XY(7»1 )*C03  004100 

XSNRI I )=100.*CO  004110 

YSNR( I )=100.*SO  004120 

XSN7(I)=(100.+SIGN7(I) >  *CO  004130 

YSN7( I )  =  (100.  +  SIGN7(  I  )  )*SO  004140 

XTNR( I )=40.*C0  004150 

YTNRI I )=40.*SO  004160 

XTN7( I )=(40.+TAUN7(I ) )*CO  004170 

YTN7( I )=<40.+TAUN7( I ) ) *SO  004180 

26  CONTINUE  004190 

C   THIS  IS  THE  END  OF  PART  8  AND  THE  START  OF  PART  9.   PART  9  CONTAINS  004200 

C   THE  STATEMENTS  WHICH  CALL  THE  SUBROUTINE  TO  PLOT  THE  STRESS  004210 

C   DISTRIBUTION  PATTERN.  004220 

CALL  DRAW(628»XSNR»YSNR»1.0»LD»ITITLE»50.»50.,4,4,2t2»8»8»0»LAST)  004  2  30 

CALL  DRAW(628»XSN7,YSN7,2t0,LB,ITITLE»50.,50.,4,4,2.2,8,8»0,LAST)  004240 

CALL  DRAW( 12 ,XSXX , YSXX ,3, 1 ,LD, I TI TLE, 50 . ,50. ,4,4, 2 , 2 ,8, 8 , 0, LAST )  004250 

ITITLE(7)=8H  TAU  VS.  004260 

CALL  DRAW(628,XTNR,YTNR,1 ,0 ,LD, I T I TLE »20 . , 20. ,4,4,2 , 2 , 8, 8 , 0 , LAST )  004270 

CALL  DRAW ( 628, XTN7 ,YTN7  ,2 ,0 ,LB , I T I TLE .20. ,20. ,4,4 , 2, 2,8,8,0, LAST)  004280 

CALL  DRAW( 12 »XTXY, YTXY  ,3 , 1  ,LD, I T I TLE ,20 . ,20. , 4,4, 2 ,2 ,8 , 8 » 0 , LAST )  004  290 

C   THIS  IS  THE  END  OF  PART  9  AND  THE  START  OF  PART  10.   PART  10  COMPUTES  004300 

C   THE  THREE  FORCES  AND  THE  THREE  MOMENTS  AT  EACH  RING  WHICH  ARE  USED  004310 

C   FOR  THE  EQUILIBRIUM  CHECK.  004320 

804  XFOR(K)=2.*3.14159*ZX(  1,1)*RA(K)*TH(K.)  004330 

YFOR(K)=-3.14159*XY( 2 , 1 ) *RA ( K ) *TH ( K )  004  340 

ZFOR (K ) =3.14 159*XY< 3,1  )  *RA ( K. ) *TH( K )  004  350 

XMOM(K)=2.*3.14159*XY( 1  ,  1 ) *RA { K > *RA ( K ) *TH< K )  004  360 

YM0M(K)=3.14159*ZX<2»1  )  *RA < K ) *RA ( K ) *TH ( K )  004  370 

ZMOM(K)=-3.14159*ZX( 3,1 ) *RA ( K ) *RA ( <) *TH ( K )  004  380 

85  CONTINUE  004390 

PRINT  5  55,(XFOR(K) ,YFOR(K> ,ZFOR(K) ,XMOM( K ) ,YMOM ( K ) ,ZMOM(K) ,<=1,8)  004400 

55  5  FORMAT ( ///2X7HXFOR < K ) 5X7HYFOR ( K) 5X7HZFOR ( K ) 5X7HXMOMI K ) 5X7HYMOM( K ) 5  004410 

1X7HZM0M(K)///(6E12.5) )  004420 

END  004430 

C   THIS  IS  THE  END  OF  PART  10  AND  THE  START  OF  PART  11.   PART  11  IS  THE  004440 

C   SUBROUTINE  GAUSS  2.  004450 
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SUBROUTINE  GAUSS2 t N,M , EP , A  ,B . X , KER )  00^10 

DIMENSION  Adl.H)  ,B(11.1)  tX(ll,l)  00020 

npm=n+m  00030 


DO  2  <  =  l.M 
I  =  N  +  K 


00040 

00050 


DO  2  J  =  1,N  00060 

2  A( J, I )  =  B( J,<)  00070 

10  DO  34  L=1,N  00080 
KP  =  0  00090 
z=0«0  00100 
DO  12  K  =  L,N  00110 
IF(Z-ABSF(A(K.L) ) ) 11 tl2.12  00120 

11  Z=ABSF< A(K,L) )  00130 
<P=K  00140 

12  CONTINUE  00150 
IF(L-KP)13,20.20  0O160 

13  DO  14  J=L,NPM  00170 
Z=A(L.J)  00180 
A(L,J)=A(KP.J)  00190 

14  A(KP,J)=Z  00200 
20  IF(ABSF(A(L,L) )-EP)50, 50,30  00210 

30  IF(L-N)31,40,40  00220 

31  LP1=L+1  00230 
DO  34  K=LP1,N  00240 
IF(A(K,L>  132,34,32  00250 

32  RATIO=A(K»L)/A<L»L>  00260 
DO  33  J=LP1,NPM  00270 

33  A(K,J)=A(K, J)-RATIO*A(L»J)  00280 

34  CONTINUE  00290 

40  DO  43  I =1 »N  00300 
I  I  =N  +  1- I  00310 
DO  43  J=1,M  00320 
JPN=J+N  00330 
S=0.0  00340 
IF( I I-N)41,43,43  00350 

41  IIP1=II+1  00360 
DO  42  K=I IP1 ,N  00370 

42  S=S+A(  I  I.K)*X(K»J)  00380 

43  X( I  I  ,J)  =  (A(  I  I  ,JPN)-S)/A( I  I ,1  I  )  00390 
KER=1  00400 
RFTURN  00410 

50  KER  =  2  00420 

END  00430 

C   THIS  IS  THE  END  OF  PART  11  AND  THE  START  OF  PART  12.   PART  12  IS  THE    004460 

C   SUBROUTINE  DRAW.  004470 

SUBROUTINE  DRAW  (NUMPTS,  X,  Y,  MODCURV,  ITYPE,  LABEL,               00000000 

1  ITITLE,  EXSCALE,  YSCALE,  IXUP,  IYRIGHT,  00000010 

2  MODEXAX,  MODEYAX,  IWIDE,  IHIGH,  IGRID,  00000020 

3  LAST)  30 
C  0  40 
C  0  50 
C  A  GENERAL  CURVE  DRAWING  AND  POINT  PLOTTING  SUBROUTINE  0000006 
C  PROGRAMMER  J.  R.  WARD  00000070 
C  DATF  FEB.  1964,  RFVISED  JUNE  1965  00000080 
C  SYSTEM  FORTRAN  60  00000090 
C          OUTPUT        LOGICAL  TAPE  NUMBER  8  00000100 
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C  NOTE          ASTERISKS  MA^K  CHANGES  FJK  EOKT^  .'•  \  63              DC00011 

C  000*  12C 

C      INPUT  ARGUMENTS ^t>"13 

C  "0""   14" 

C  1.  NU'-iP1"0,      .\."<DF-  OF  JOINTS  TO  ^c  PLOTTED.  THIS  M'JST  ALW  A v  S  ?  0° ?  0  1  50 

r  BF.  AT  LEAST  ?,  ANO  MUST  NOT  EXCEED  30  rQR  JOINT  00^0160 

C  PLOTTING*  2»  900  cjR    CURVE  DRAWING.              00000170 

C  CCO*  180 

C  2.  X            ARRAY  Of  X-OROINaTES.  DIMENSION  AT  LEAST  EQUAL   00000190 

C  TO  NUMPTS  AND  NOT  MCRF  THAN  >'""v  [N  CALLING       O0C0O20O 

C  PR'JGKAv.                                           00000210 

C  COC*  220 

C  3.  Y           ARRAY  OF  Y -CRD  I  NATES.  DIMENSION  AS  FOR  THE       O0OC0230 

C  X  ARRAY  IN  THF  CALLING  PROGRAM.                    00000240 

C  000   250 

C  4.  MODCURV     CONTROLS  THF  NUMBER  OF  CURVES*  AND/OP  SETS  0e          00000260 

C  POINTS*  ON  FACH  GRAPH.  THF  COOES  ARE              CC~Cn2?0 

C  0     ONLY  ONF  PLOT  ON  THIS  GRA"H           0C0OO280 

C  1     FIRST  PLOT  ON  MULTI-PLOT  GRAPH        00^00290 

C  2     INTERMEDIATE  PLOT  ON  MULTI-PLOT  GRAPHOO000 300 

C  3     LAST  PLOT  ON  MULTI-PLOT  GRAPH.        00000310 

C  000   320 

C  5.  ITYPE       CONTROLS  THE  TYPE  Or  PLOT.  THE  CODES  ARE         00000330 

C  0     STRAIGHT  LINES  JOIN  SUCCESSIVE  PO I  NT SOOOC 0 340 

C  (STANDARD  CURVE  DRAWING)               00000350 

C  1     POINTS  PLOTTED  WITH  CROSS  (X)         OOOO036O 

C  2     POINTS  PLOTTED  WITH  PLUS  (+)          00000370 

C  3     POINTS  PLOTTED  WITH  SQUARE            00000380 

C  A     POINTS  PLOTTFO  WITH  DIAMOND            O0OOQ39O 

C  5     POINTS  PLOTTED  WITH  TRIANGLE          00000400 

C  WHFN  POINTS  ARE  BEING  PLOTTFD  !ITYPF=1  THRU.  5),000OO^io 

C  THE  POINTS  ARE  NOT  CONNECTED.                     00000420 

C  000   430 

C  6.  LABEL       IF  A  CURVE  IS  BEING  DRAWN  ( ITYPE  =  0).  LABEL  IS  00000440 

C  A  4-CHARACTER  BCD  CURVE  IDENTIFIER  WHICH  WILL  6E00000450 

C  REPRODUCED  AT  THE  END  OF  THE  CURVE.  FOR  EX  AMPLE , 00000460 

C  LABEL  =  4H  ONE  .  IF  POINTS  ARE  HEING  PLOTTED,    O0000470 

C  LABEL  IS  AN  8-CHARACTFR  IDENTIFIER.  THE  FIRST  4  ^0000480 

C  CHARACTERS  ARE  REPRODUCED  WITH  THF  FIRST  PLOT TED00000490 

C  POINT,  AND  THE  LAST  4  WITH  THE  LAST  POINT.  SET   00000500 

C  TO  BLANK  ANY  UNWANTED  CHARACTERS.                 OO000510 

C  000-  520 

C  7.  ITITLE      AN  ARRAY  OF  TWELVE  8-CHARACTER  BCD  WORDS.        O0000530 

C  THE  FIRST  SIX  WORDS  WILL  BE  REPRODUCED  AS  THE    00000540 

C  FIRST  LINE  OF  GRAPH  TITLE,  AND  THF  LAST  SIX      00000550 

C  WORDS  WILL  FORM  THE  SECOND  LINE.  TH"7  TITLE  MUST  0000C560 

C  INCLUDE  THE  USERS  JOB  IDENTIFICATION.  DIMENSION  00000570 

C  12  IN  CALLING  PROGRAM*  AND  SET  TO  liLANK.  ALL      00000580 

C  UNWANTED  CHARACTERS.                              00000590 

C  000   600 

C  8.  EXSCALE     X-SCALE  (UNITS/INCH)  AS  POSITIVE  FLOATING  POINT  00000610 

C  VARIABLE  WITH  ONE  FIGURE  SIGNIFICANCE.  SET  TO    00000620 

C  ZERO  FOR  AUTO-SCALE.                              00000630 

c  000   640 

C  9.  YSCALE      Y-SCALE  (UNITS/INCH)  AS  POSITIVE  FLOATING  POINT  00000650 

r  VARIABLE  WITH  ONE  FIGURE  SIGNIFICANCE.  SET  TO    O0000660 
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C  ZERO  FOR  AUTO-SCALE.                              "000067" 

C  COO   680 

C  10.  IXUP        X-AXIS  OFFSET  FROM  BOTTOM  OF  G^APH  I  iM  INCHES.    00000690 

C  THIS  MUST  NOT  EXCEED  I H I GH ,  AND  MUST  NOT  BE      00000700 

C  NEGATIVE.                                          00000710 

C  COO*  720 

C  11.  IYRIGHT     Y-AXIS  OFFSET  FROM  LEFT  OF  GRAPH  IN  INCHES.  THIS00000730 

C  MUST  NOT  EXCFED  IWIDE.  AND  MUST  NOT  BE  NEGAT I VE .00000740 

<~  000   750 

C  12.  MODEXAX     MODE  OF  X-AXIS  OFFSET.  SEE  CODES  BELOW.          0000076" 

C  COO   770 

C  13.  MODEYAX     MODE  OF  Y-AXIS  OFFSET.  THE  CODES  ARE  AS  FOLLOWS  00"0"780 

C  0     COMPUTED  OFFSET.  HOLDING  ORIGIN  ON    00000790 

C  GRAPH.  THE  CORRESPUNDING  IXUP  OR      00000800 

C  IYRIGHT  IS  IGNORED                     00000810 

C  1     COMPUTED  OFFSET,  WITH  ORIGIN  OFF  THE  00000820 

C  GRAPH  IF  APPROPRIATE.  THE  CORRESPOND-00000 830 

C  ING  IXUP  OR  IYRIGHT  IS  IGNORED.  USE   00000840 

C  ONLY  WITH  AUTO-SCALE                   00000850 

C  2     AXIS  OFFSET  AS  SPECIFIED  BY  IXUP  OR   00000860 

C  IYRIGHT.                               00000870 

C  000*  880 

C  14.  IWIDF       WIDTH  OF  GRAPH  IN  INCHES.  THIS  MUST  NOT  EXCEFD   00000890 

C  NINE.   ZERO  WILL  BE  READ  AS  EIGHT  INCHES.        00000900 

C  000   910 

C  15.  IHIGH       HEIGHT  OF  GRAPH  IN  INCHES.  THIS  MUST  NOT  EXCEED  00000920 

C  CIFTEEN.  ZERO  WILL  BE  READ  AS  EIGHT  INCHES.      00000930 

C  000   940 

C  16.  IGRID       IF  SET  TO  1,  A  ONE  INCH  BY  ONE  INCH  GRID  WILL    00000950 

C  BE  SUPERIMPOSED  ON  THE  GRAPH.                     0000"960 

C  0000  970 

C  17.  LAST        INDICATES  TO  CALLING  PROGRAM  WHETHFR  LAST  PLOT   "0000980 

C  WAS  COMPLETED  SUCCESSFULLY.  THE  CODES  ARE        0000"990 

C  0     LAST  PLOT  COMPLETED  SUCCESSFULLY      000n1000 

C  1     LAST  PLOT  NOT  SUCCESSFUL               00001010 

C  2     LAST  PLOT  NOT  SUCCESSFUL  AND  NO       00001020 

C  FURTHER  GRAPH  OUTPUT  WILL  BF  A TTEMPT-00001 030 

C  ED  UNTIL  MODCURV  IS  NEXT  ONE  OR  ZERO  00001040 

C  3     DRAW  WAS  ENTERED  WITH  MODCURV  NOT     00001050 

C  EQUAL  TO  ONE  OR  ZERO  WHILE  THE  ERROR  00001060 

C  LOCK-OUT  WAS  SET.                      O0001C70 

C  THIS  ARGUMENT  MUST  ALWAYS  BE  A  NAME  IN  THE  CALL  00001080 

C  STATEMENT.  NEVER  A  NUMBER.                        00001090 

C  00001100 

C  00001110 

c      N0TF  00001120 

C  ALL  ARGUMENTS  FROM  NUMBER  7  THRU.  NUMBER  16  ARE  IGNORED  WHEN   00001130 

C  MODCURV  IS  FITHER  2  OR  3.  HOWEVER,  ARGUMENTS  MUST  NEVER  BE     00001140 

C  OMITTED  FROM  THE  CALLING  STATEMENT.  IT  IS  MERELY  THEIR  VALUES  00001150 

C  WHICH  ARE  THEN  IRRELEVANT.  ARGUMENTS  MAY  Bt  LISTED  BY  NAME  OR  00001160 

C  VALUE  IN  THE  CALL  STATEMENT.  NO  VALUE  IN  THE  CALLING  PROGRAM   00001170 

C  WILL  BE  ALTERED  BY  THIS  SUBROUTINE.                             00001180 

c  00001190 

C      RFFFRFNCF  —  ^"001200 

C  THE  BINARY  TAPE  FORMAT  REQUIRED  BY  THE  OFF-LINE  PLOTTER  IS     00001210 

C  DFSCPIBFD  IN  THE  WRITFUP  OF  THE  CDC  160A  GRAPH  PLOT  PROGRAM    00001220 


C  n  p  '.' 


IN- 


JECT)      Q  V       TUJC 

3 L  A T I  ON    •'  -  a 


L  E  <  1  i  >  ,  J  X  T  I  T ( 1 2  )  i  J  Y  T  I  T  (  1 2  1  » 
,  ICURVI460).  JGMD125)*  ICvV 


-3. 


007^ 


9071 


DIMENSION  X(9-'0),  H900).  HI 
L  T  !  T  L  F  (  i  <♦  )  »  <  A  K  15(5 
JJT ITLEC 121 
JPOINT  =  ITYPC 

CON(  ICONTRL  =  4':00oc). 

CON(ICURV3  =  2777377"'20?'>202On  ,  rc'JRV*  =  p  1 n  '■ 
REPLACE  'a!th  niT<  STATEyENT  !N  FO'TPAM  6; 
PUT  ITC<iT  =  *">  IN  TAT".  STAT^m^NT. 
IF  (JTEC.T  -  735t6912)   9070 ,  9n  7j  ,9  r  7n 
ITFST  =  0 
JTEST  =  73546912 
CONTINUE 

REMOVE  ABOVE  NONSENSE  IN  FORTRAN  63. 
CHECK  PREVIOUS  OPERATION  OF  ROUTINE.  IF  ANY.  CODES  ARE 
IF  PREVIOUS  GRAPH,  IF  ANY,  COMPLETED 
IF  PREVIOUS  GRAPH  NOT  COMPLETED 


ITEST 
ITEST 
ITEST 


1001 
1002 

1003 
1004 


1005 
1006 
1007 
1100 


IF  ERROR  FOUND  WHILE  MODCURV  WAS  CNF.  ,  OR  IF 
^OOCUPV  WAS  ILLEGAL. 

IFUTEST  -   2)looofinni,ioon 

IF  (MODC'IRV)  lAr3,10o;>,ioo? 

ITCST  =  0 

GO  TO  1000 

IFIMODCURV  -  1)1004.1002.1004 

LAST  =  3 

RETURN 

SET  UP  ERROR  RETURN  ROUTINE.  ENTRY  AT  STATEMENT  1005. 
IF( I  TEST) 1009.1 006 » 1009 
IF(MODCURV) 100  7,1003,1007 
PRINT  1100 

NO  FURTHER  GRAPH  OUTPUT  UNTIL  MOOCJPV  NFXT  IS  ZERO 


1008 
11C1 


1009 
1010 
1011 


100" 
1 

ion 

2 

9000 
9100 

PI!!! 


2) 1010,  100P.1010 
3  )  1 0  0  7  ,  1 0 1 J  »  1  C  C  7 


FORMAT  (  59H 

ONE.  ♦/) 
ITEST  =  2 
LAST  =  2 
RETURN 
PRINT  1101 

FORMAT  (  30H  THIS  PLOT  WILL  NOT  6E  OUTPUT.  ,/) 
LAST  =  1 
RETURN 
IFIMODCURV 
IF(MODCURV 
ITEST  =  0 
GO  TO  1008 

CHECK  LEGALITY  nt 
IFtNUMPTS  -  2)1,2,? 
PRINT  100 

FORMAT  (/,  32H  NUMPTS  MUST  N0T  BF  LCSS  THAN  2, 
GO  TO  100  5 

I<~  (  IPOINT  )9flrn,q^o^,9001 
PRINT  9100 

FORMAT  (/,  15H  ILLEGAL  [TYPE.  ) 
GO  TO  10C5 
IF(  IPOINT  -  5)9002.9002.9000 


INPUT  ARGUMENTS. 


3CC 01250 

"::->i26? 

~000127~ 
" ~  0O123O 
0C0C1290 
3 0  0  0 1 300 
00001313 
^  0  ""M  3  ?  A 

0000  1  330 

00001  3^0 
nnT>\  3^0 

0000 1 360 
00001 370 
00001 360 
00001390 

ooooi^oo 

00001410 
OCC01420 
00C01430 
O00P144Q 
OOOO 1450 
00001460 
00001470 
O0001480 
00001490 
00001500 
00001510 
00001520 
000^1530 
00001540 
00001550 
OR00001560 
000C1570 
00001580 
00001590 
00001600 
00001610 
00001620 
00001630 
00001640 
00001650 
00001660 
00001670 
0000168'*' 
00001690 
00001700 
00001710 
00001720 
O00C173O 
00001740 
00001750 
00001760 
00001770 
00001780 
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0002  IF(N'.JWPTS  -  30)3,3,9^03  00001790 

ono^  OD]f,<r    qini  OOnnisno 

0  1 01  FORMAT  (/,  '46H  NUMPT*  MUST  NOT  CXC!!E''  ^0  FOR  °G'I\T  PLOTTING.  |     onn^isiO 

GO  TO  10r5  ^0001820 

qO04  [F(Nl.'MPT5  -  ann  )■},■»,  ono*  ^nr,1830 

9005  PRINT  9  102  000018^0 

9102  FORMAT  (/.  2SH  NUMPTS  MUST  NOT  EXCEED  910,  )                        00001850 

GO  TO  1005  00001860 

3  IX  =  1HX  00001870 
IY  =  1HY  00001880 
AMAXX  =  -0.2F+10O  00001890 
AMAXY  =  -0.2E+100  00001900 
AMINX  =  +O.2E+100  00001910 
AM  I  NY  =  +o.?E+100  OC001920 
00  1020  1=  1, NUMPTS  00001930 
AMAXX  =  MAX1FIXI  I  )  ,AVAXY)  00001940 
AMAXY  =  MAX1FI  Y(  I  )  , AMAXY)  00001950 
Aminx  =  MIN1F<X( I ) tAMINX)  00001960 

1020  AMINY  =  MIN1FI  Y(  I  ) ,AMINY)  00001970 
AMAXA  =  MAX1FI  A3SFI  AMAXX)  »  AdSF(AMAXY),  A3SF ( AM  I  NX ) ,  ASSF(AMINY))  00001980 
IF(AMAXA  -  1.0E+99) 1  022 » 1022 , 1 A2 1  00001990 

1021  PRINT  1102  00002000 

1102  FORMAT  (/,  58H  NO  X  OR  Y  VAL'JF  MAY  EXCEED  1.0F  +  99  IN  ABSOLUTE  MAGN00002010 
1ITUDE.  )  00C02020 

GO  TO  100*;  00002030 

1022  IF(ABSF( AMAXX  -  AMINX)  -  1  .  OE-97 ) 1 023 , 1 025 , 1 02 5  OQ002040 

1023  IF(ABSF( AMAXY  -  AMINY)  -  1 .OE-Q7 ) 1 024 , 1 025 » 1 02 5  00002050 

1024  PRINT  1103  00002060 

1103  FORMAT  (/,  38H  ALL  POINTS  HAVF  THE  SAME  COORDINATES.  )  O0O02070 
GO  TO  1005  00002080 

1025  IF(  ITEST)4,7,4  00002090 

4  IFIMODCURV  -  2)5,240,5  00002100 

5  IFIMODCURV  -  3)6*240.6  00002110 

6  PRINT  101  00002120 

101  FORMAT  </,  17H  ILLEGAL  MODCURV.  )  00002130 
GO  TO  1005  00002140 

7  Ic (MODCURV) 6 ,9,8  00002  150 

8  IFIMODCURV  -  1)6»9,6  00002160 

9  IF(  IWIDF)  10,1 ]  ,1 2  00002170 
in  ITIT  =  5HIWIDE  00002180 

PRINT  102,  ITIT,  ITIT  00002190 

102  FORMAT  (/,   9H  ILLFGAL  ,A5,29H.  GRAPH  WILL  3E  PLOTTED  WITH  ,A5,    00002200 
1              5H  =  8.  ,/)  00002210 

11  JWIDF  =  8  00002220 
GO  TO  14  00002230 

12  IFUWIDE  -  9)13,13,10  00002240 

13  JWIDE  =  IWIDE  00002250 

14  IF( IHIGH) 15,16, 17  00002260 

15  ITIT  =  5HIHIGH  A0002270 
PRINT  102,  ITIT,  ITIT  00O02280 

16  JHIGH  =  P  00002290 
GO  TO  19  0CO02300 

17  IFMHIGH  -  15)18,18,15  00002310 

18  JHIGH  =  IHIGH  00002  320 

19  NODFXAX  =  MODEXAX  00002330 
IF(MODEXAX)20,27,21  00002340 
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2'-     I  T  I  T  =  8HM0DFXAX,  "D^C^SSO 

PRINT  104,  I  T  I  T  •  IX  00102361 

104  FORMAT  </,   9H  ILLEGAL  »A8*  32w  GRAPH  a  ILL  aS  PLOTTED  WITH  MODE*   000O237C 

1          Alt   7HAX  =  •"  ,  ,/)  0OOO2!3^ 

NODEXAX  =  n  0000239C 

r^O  to  27  Ov(",0240" 

?!  IF(MODEXAX  -  2)27t22»2C  00002410 

22  !F(IXMP  -  JHIGM)  24,24,2'.  00002420 

23  ITIT  =  8HIXI.IP.  00002430 
DOINT  lr4,  inT»  U  ~0002440 
NODFXAX  =  1  00002450 
GO  TO  27  0000246C 

24  IF  (  IX'JP)23,26.26  10002470 

26  JX'JP  =  IXUP  00002480 

27  NODEYAX  s  MODEYAX  00002490 
IF(MODEYAX)28i35»29  0000250O 

2«  ITIT  =  9HMODcYAX.  "0C02510 

dqint  1^4,  ITIT,  IY  0C0P2520 

NjDOFYAX  =  0  00002530 

GO  TO  35  00002540 

20  IF(M.ODEYAX  -  .")35,30,28  00002550 

31  irfjYPIGHT  -  JWIor )32.32.31  10102560 

31  ITIT  =  8HIYRIGHT,  00002570 
PRINT  104,  ITIT.  IY  00002580 
NODEYAX  =  C  00002590 
GO  TO  35  00002600 

32  IFUYRTGHT)  31.34,34  00002610 

34  JYPIGHT  =  IYRIGHT  00002620 
C  INITIALIZE  PRIOR  TO  SCALING  AND  AXIS  LOCATING.  00002630 
C           IFLAG  =  0  FOR  PASS  WITH  XHATA.  IFLAG  =  1  FOR  PASS  WITH  YDATA.  00002640 

35  DO  2235  IOTA=l»12  00002650 
2235  JJTITLEt IOTA)  =  ITITLE(IOTA)  10002660 

IFLAG  =  0  00002670 

B^TA  =  0.  00002680 

SCALE  =  EXSCALE  00002690 

IAXIS  =  JYRIGHT  10002700 

MODE  =  NODEYAX  00002710 

ISIZE  =  JWIDE  00002720 

IXY  =  IX  00002730 

IYX  =  IY  01002740 

A^AX  =  AMAXX  00002750 

AMIN  =  AMINX  00002760 

GO  TO  52  00012771 

5i  IFLAG  =  1  00002780 

RFTA  =  0.  00002790 

SCALE  =  YSCALE  00002800 

IAXIS  =  JXUP  00002810 

MODE   =  NODEXAX  00002820 

ISIZE  =  JHIGH  00002830 

AMAX  =  AMAXY  00002840 

AMIN  =  AMINY  00002850 

TXY  =  IY  00002860 

IYX   =  IX  00002870 

C          CHECK  SCALE  AND  GO  TO  F1XFO  OR  AUTO  SCALE  ROUTINES.  00002880 

5?  IF(SCALF)53,59,56  00002890 

53  PRINT  114,  IXY,  IXY  00002900 

59 


14 


K  7 


10  3' 


10 

11 

1^ 
10 
10 
10 

10 

10 

10 
10 


31 
04 


32 

3  3 
3  0 


36 
37 

38 
40 


5  ° 
60 

64 

6? 

I  16 


6B 


69 

7o 

71 

.  18 

74 
75 
76 

77 


78 
70 


FORMAT     (/,       9-i 

T^     ,A1,     7H-crAL 

",  r\     t  g     =,  q 
ryD^ccc     cjyri-     c 
rF(SrA|_-     -     ],rc 
Ic (  S^*LC     -     1  ."r 
CALL    r  r  '■  L  r  !  T  (  ?  ' 

5C4^  =  CjrT-:!3< 
fucri;  A  N  0  r 
C  A  5  ~  .     !  T  A  C-. 

I  c  (  MOOF    -     1  )  1  'i? 

!  T  A  G    =  0 

GO    TO    203 

PR  I  NT     1 lr'4    ,     l  y 

FORMAT  (/,  5H  v 
(  VJTO-SCALE)  . 

r.r\     t  Q     c  Q 

I  F  (  A  P  "  c  (  A  V  A  X      - 

I  F  (  A.RSF  (  ftVAX  )      - 

IF<  Av-.xi  10  3  6,1" 
TAXIS  =  I  <-.  I  Z  E  /  3 
GO  TO  103  0 
I  AXIS  =  ISIZE 
GO  TO  1030 
I  AX  IS  =  0 
GO  TO  1030 
IF(SIGNF( 1. ,AMA 
ASI7E  =  ISIZE 
I  AX  IS  =  -AMIN/I 
GO  TO  103O 

AUTO  SCALE 
IF(M0OF  -  1)60, 
AMAX  a  MAX1F(0. 
AM  IN  =  MIN1F(0. 
IF(AP.SF(AMAX  - 
PRINT  116,  IXY, 
FORMAT  (/,   5H 
IF  THF  ,A1»29H 
0,0  TO  100  5 
ASIZE  =  ISIZF 
SCALF  =  (AMAX  - 
iSO  TO  8  3 
IF(APSF(AMAX  - 
IF ( ARSF ( AMAX )  - 
PRINT  118,  IXY 
FORMAT  (/»   5H 
GO  TO  1005 
IF (AMAX  -  l.OE- 
IF( ISIZE  -  IAXI 
SCALE1  =  0. 
CO  TO  7« 
AXIS  -  I  A  X  I S 

A  S  I Z  F  =    I «;  I  Z  E 
SCALE]  =  AMAX/I 
[  F  (  A  M I  N  +  I  •  0  F  - 

IF(   T  AX  t  c  )  P.  1  ,80, 


!LL5G^L 
r  .     .  /  ) 


.  1  ^  u  c  r  a  1 


•      :  ■-! 


ALE.I  '-"*_■      »rA^T-»3,i  1 
l-"1.**  ISC'-L  1  " 

rjypijT  -     AXIS    LC'"  A  T  I  ON     '  c     Nc  0~  S  S  A'^  v  .     c  '  x  -  ~*     c  ^  a  ' 

=     "'     ;~     0'-' !  G  I  ■*      "\     '"A~u     "' •?     1      i  c     iT     IS     vj^rjctr 
7,inii  .  !  o  -j  •• 


>  ->  c  -J  ■> 

'206^ 

>  2  <3 "»  0 

I  jj  r  Q  1 

1  ?  ::  -  n 

■  i  ->  a-, 


?  0  0  0 


X  ,      IXY,      r  X  Y 

0DE,A1,24HAX     -"JST    NOT    3E     1    UNLFS5     »  A  W57HS.CA! 

GRAPH    'vILL     -ic     PLO^TFO    WITH    AUTO    ,A1,     7H-SCAL! 

f\MTty)     -    l.Qc-97)  1033, 1038,  1^38 

I/)--5)?)  10  34,1*39,103° 
T  4  ,  1  )  3  7 


3030 


X!     -    SIGNF< 1.,AMIN) U040, 1039,1340 
AMAX    -    AMIN)*A$IZE    +0.5 

ROUTINE. 

64,69 

,  AMAX) 

,  A  m  I  N  ) 

AMIN)  -  1.0F-97)65»68,68 

IXY,   IYX 
ALL  , A1,47H  VALUES  EQUAL.  AUTO  SCALE  POSSIBLE 
VALUES  ARE  NON-ZFRO  AND  Y0DF.A1,  7HAX  =  2.  ) 


AMINI/ASIZE 

4.MIN)     -     1.0E-97)70»74,74 

1.0F-97)  71 ,74, ?4 

ALL  ,A1,38H  VALUES  ZERO.  AUTO  SCALE  NOT  POSSI! 

97)  76,75,75 
5)77,76,77 


A  S ! / F  -  AXIS) 
07) 79,79,80 

81 


*> "~  *~\  •>  -5  ,^  -7  >\ 

0 .0  0  0  3  ^  p  n 

000/)3  09^) 

0C^C3100 

00 003110 

0  00  0  3120 

00003  130 

'"•0003140 

0  0  0  0  3150 

0  0  00  3  16  0 

00003 1 70 

CO 0  03180 

00003190 

r>  0  0  0  3  2  0  0 

^^^0"3  77^ 

00003  220 

00003230 

00003240 

00003250 

00003260 

ONLY 

0  000  3  2  70 

00003 280 

00003290 

00003300 

00003310 

ppom  3  -pr\ 

00003330 

o00033''O 

00003350 

LE. 

100003360 

00003370 

00003380 

00003390 

00003400 

00003410 

00003420 

00003430 

0000^440 

O0003450 

00003460 

60 


8^  SCALC2  =  '".  0"°034?o 

GO  TO  82  00C03490 

81  AXIS  =  I  AXIS  00003490 
SCALE2  =  -AMIN/AXIS  ?C0O35C0 

82  IFISCALF1  +  SCALC2  )  1  994  ,  1932  »  1  9<?4  00003510 
108?  PRINT  1983,  IYX.  IYX  "0003520 

1983  FORMAT  </»  56H  NONE  0*  THE  PLOT  LIES  O'l  THE  GRAPH  WITH  THIS  SPEC ! F000035 30 
1IFD  »A1»47H-AXIS  LOCATION.  SRAPH  WILL  RE  PLOTTED  WITH  ••'C0E»A1,     "0^03540 

?  7hax  =  n,      ,/>  oonn35co 

MODE  =  0  00003560 

GO  TO  60  00003570 

1984  SCALF  =  MAX1FI SC*LF1,  SCALF2)  00003580 

83  CALL  SCM-FITtSCALF,  I  SC  AL  1  0  ,  F  -\  rTCR  ,  3  )  00003c90 
IF(FACTOR  -  5.05! 35,85,84  00003600 

84  FACTOR  =  1  0C003610 
ISCAL10  =  ISCAL10  +  1  00003620 
GO  TO  90  00003630 

85  IFtFACTOR  -  2.02)87,87,86  00003640 

86  FACTOR  =  5  CC003650 
GO  TO  90  om~?6ic)~ 

87  IF(FACTOR  -  1.01)89,89,88  00003670 

88  FACTOR  =  ?  00003680 
GO  TO  90  00003690 

89  FACTOR  =  1  00003700 

90  SCALF  =  FACTOR*10.**TSCAL10  00003710 
C          COMPUTE  AXIS  LOCATION  IF  NECESSARY.  AUTO  SCALE  CASE.  00003720 

IFtMODE  -  1)92,91,93  00003730 

91  IF(SIGNF( 1.  ,AMAX)  -  S 1 GNF ! 1 . , AMI N )  ) 92 , 94  ,  92  00003740 

92  I  AXIS  =  -AMIN/SCALE  +  0.5  00003750 

93  ITAG  =  0  00003760 
GO  TO  203  00003770 

94  IF( AMAX)95,95.2O0  00003780 

95  IAXIS  =  ISIZE  00003790 
BFTA  =  -AMAX/SCALF  0000380O 
!F<eETA  -  1. £+12)99,99, 96  00003810 

96  PRINT  1?0,  IXY  00003820 
12C  FORMAT  (/,  15H  THE  OR/GIN  OF  ,A1,  43H  CANNOT  BE  OFFSET  MORE  THAN  100003830 

l.OF+12  INCHES.  )  00003840 

GO  TO  1005  00003850 

99  IRETA  =  BETA  +  0.5  00003860 

RFTA  =  -IBFTA  00003870 

C  BETA  IS  THE  r-i'JMBPR  OF  INCHES  OF  ORIGIN  SUPPRESS  I  ON, POS I T I VE  IF00003880 

C           TRUC  ORIGIN  1^  RFLOW  OP  TO  LFFT  OF  THC  GRaPh.  00003890 

IF(BETA  +  1.197,97,93  00003900 

97  ITAG  =  1  00003910 
GO  TO  ?n?  00003920 

200  IAXIS  =  0  00003930 
BFTA  =  AMIN/SCALF  00003940 
IFIBFTA  -  l.E+12)201 ,201,96  00003950 

201  1BETA  =  BFTA  +  0.5  00003960 
BETA  =  IPFTA  00003970 
IFfBETA  -  1.  193,202,202  00003980 

202  ITAG  =  1  00003990 
C  RELEASE  RESULTS  TO  REMAINING  PART  OF  PROGRAM.  START  SECOND  00004000 
r           PASS  FOR  Y  V'ALIJFS  IF  NOT  YFT  COMPUTFD.  00004010 

203  IF( I FLAG) 205* 20 4* 2 05  00004020 

61 


214 


20  5 


206 


216 


c^Lpy 

IX- ATT 

r  x  s  c  l  o 

IXAXIS 

ITAGX 

ISIZEX 

PETAX  = 


=  <^CALE 

=  FACTOR 
=  ISCALin 
=  IAXIS 
=  ITAC- 
=  I  SIZE 
BFTA 


9206 

207 
208 

20Q 

2in 


211 

212 
213 

214 
?15 


GO  TO  50 
^TAY  =  BFTA 
SCALFY  =  SCALE 
IYFACT  =   FACTO'' 
I  v$r jn  a  t  crsL  1  n 
IYAXIS  =  IAXIS 

NOW  WRITE  RECORDS. 
ITAGY   =  I  TAG 
ISIZEY  =  ISIZE 

THIS  COMPLETES  CALCULATI 
RECORDS.  FIRST,  THE 

=  em     x- 

=  8HSCALF  = 

=  ICODc< SCAL^X) 

=  8H  UN  ITS/ I 

=  BHNCH. 

=  8H1      Y- 

=  PHSCAL11  = 

=  rCODE(SCALFY) 

=  8H  UNITS/ I 

=  8HNCH. 
1=6.11 

=  8H 

=  8H 

IF( ITAGX)211,211,207 
IF(BETAX)208»208.209 
JXTIT( -,  =  8H    ADD  - 
GO  TO  210 

JXTIT( 7)  =  8H    ADD  + 
JXTIT(8)  =  ICODP  (9ETAX»SCALFX) 
JXTIT(9)  =  8H  UNITS  T 
JXTITI 10) =  8HO  ALL  X 
JXTIT(11)=  8HVALUES. 
IF(  ITAGY)216.216»212 
IF(BETAY)213,213.214 
JYTIT(7)  =  8H    ADD  - 
GO  TO  215 


TAPE 
JXTIT(  1  ) 
JXTITf  2) 
JXTIT( 3  ) 
JXTITI4) 
JXTIT( 5  ) 
JYTIT( 1 ) 
JYTITI 2) 
JYTIT( 3) 
JYTIT(4) 
JYTIT(5) 
DC  9206 
JXTIT(  I  ) 
JYTITt  I  ) 


>N  OF    SCALE  FACTORS  ETC.  NOV 
SCALE  FACTOR  TITLES. 


JYTIT(7)  = 
JYTITI8)  = 
JYTIT(9)  = 
JYTITf 10) =8HO  ALL  Y 
JYTITI 11 )=  8HVALUFS. 
irONTt 1 )  =  ICONTRL  + 


8H    ADD  + 

ICODE  (BETAY*SCALEY) 

8H  UNITS  T 


INSERT  TITLF  SIZF  (02B)  AHPAD 
CALL  ISHIFT6  (ITITLE,  LTI1LF) 
TEST  FOR  ALL  BLANK  TITLES. 
ICHECK  =  8H 
DO  9075  1=1,6 
IF (  I  TITLE (  I  ) -ICHECK)  9074,90  75,9074 


OF  MAIN  TITLE  RECORD. 


?  *  i 

ceo 

00n 
0  CO 

ceo 
ceo 
c  re 
coo 

00*"! 

000 

coo 
coo 
coo 
coo 
coo 

GENERATE  000 

COO 
000 
000 
0  00 
000 

oco 

000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 


04030 
14141 

C4C50 
0^060 
34070 

C4C80 
04  090 
OMOO 
04  11C 
04120 
04130 
04141 
04150 
04161 
04170 
04181 
04190 
04200 
04210 
04220 
04230 
04240 
04250 
04260 
04270 
04280 
04290 
04300 
04310 
04320 
04330 
04340 
04350 
04360 
04370 
04380 
04390 
04400 
04410 
04420 
04430 
04440 
04450 
04460 
04470 
04480 
04490 
04500 
04510 
04520 
04530 
04540 
04550 
04560 
04570 
04580 
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9074  I  F ( ITT  TLE ( I  1  )  9080.9075,9^80 

9075  CONTINUE 
IT1  =  1 

ICONT(I)  =  ICONT(l)  -  1 
GO  TO  9081 

9080  IT1  =  0 

9081  DO  9085  1=7,12 

IF  (ITITLE(I)  -  ICHECK)  9084 t 9085 , 9084 

9084  IF  (  ITITLE( I  )  ) 9090 »9085 #9090 

9085  CONTINUE 
IT2  =  1 

ICONT(l)  =  ICONT(l)  -  1 
GO  TO  9091 

9090  IT2  =  0 

:  NOW  GFNERATF  AXES  RECORDS, 

9091  LFTMGN  =  0*100 
IBOTMGN  =  1*100 
I'-i  -  LFTfGN 

jK  =  IBOTMGN  +IYAXIS*100 
LH  =  ISIZFX*100 

IHL  =  LFTMGN  +  ISIZEX*100  -  107 
KAXIS(l)  =  IPACK12(IH,JH,LH,IHL> 
JHL  =  JH  -  13 
IHL2  =  -100 

IVH   =  (ISIZEX  -IXAXIS  -  1)*IXFACT 
IVH2  =  -IXFACT 

KAXISJ2)  =  IPACK12( JHL . I HL2 . I VH, I VH2 ) 
NH  =  ISIZEX 
ISH  =  8H       14 
IV  =  LFTMGN  +  IXAXIS*100 
JV  =  IBOTMGN 

KAXISO)  =  IPACK12<NH,ISH»IV»JV) 
LV  =  ISIZEY*100 
IVL  =  IV  -  3 

JVL  =  IBOTMGN  +  ISIZEY*100  -  107 
JVL2  =  -100 

KAXIS<4)  =  IPACK12(LV»IVL»JVL,JVL2) 
IVV  =  (ISIZEY  -  IYAXIS  -  1)*IYFACT 
IVV2  =  -IYFACT 
INV  =  ISIZEY 
ISV  =  8H       11 

KAXIS(5)  =  IPACK12UVV.IVV2.INV.ISV) 
I  NOW  GENERATE  CURVES. 

SCX  =  100./SCALEX 
SCY  =  100./SCALEY 
EXAXIS=  IXAXIS*100 
YAXIS  =  IYAXIS*100 
ALFTMGN  =  LFTMGN 
BOTMGN  =  IBOTMGN 

SHIFTX  =EXAXIS  -  BFTAX*100.  +  ALFTMGN 
SHIFTY  =  YAXIS  -  BFTAY*100.  +  ROTMGN 
EXSIZF=  ISIZEX*100  +  LFTMGN  +  60 
SIZEX  =  LFTMGN  -  60 
YSIZE  =  ISIZEY*100  +  IBOTMGN  +  70 
SIZEY  =  IBOTMGN  -  70 
ICURV(l)  =  0 


O00045°0 
C00046CC 
00004610 
00004620 
00004630 
00004640 
00004650 
00004660 
00004670 
00004680 
00004690 
00004700 
00004710 
00004720 
00004730 
00004740 
00004750 
00004760 
00004770 
00004780 
00004790 
00004800 
0C004810 
00004820 
00004830 
00004840 
00004850 
00004860 
00004870 
00004880 
000C4890 
00004900 
00004910 
00004920 
00004930 
00004940 
00004950 
00004960 
00004970 
00004980 
00004990 
00005000 
00005010 
00005020 
00005030 
00005040 
00005050 
00005060 
00005070 
00005080 
00005090 
00005100 
00005110 
00005120 
00005130 
00005140 
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240  IF( IPOINT 19010.9007,9010  00005150 
9007  IF(XMODF(NUMPTS»2) ) 9700 .9701 » 9700  00005160 

9700  ISWITCH  =  1  00005170 
GO  TO  242  00005180 

9701  ISWITCH  =  2  00005190 
242  INUM  =  (NUVPTS  +  1)/?  00005200 

DO  244  I=1,INUM  00005210 

CI  =  X(2*I-1)*SCX  +  SHIFTX  00005220 

C2  =  Y(2*I-1)*SCY  +  SHIFTY  00005230 

IF<  I-INUM)241, 9241, 241  00005240 

9241  GO  TO  (9242,241) .ISWITCH  00005250 

9242  C3  =  CI  00005260 
C4  =  C2  00005270 
GO  TO  9243  00005280 

241  C3  =  X(2*I)*SCX  +  SHIFTX  00005290 
C4  =  Y(2*I)*SCY  +  SHIFTY  00005300 

9243  CI  =  MIN1F(C1,EXSIZE)  00005310 
IC1=  MAX1FIC1,  SIZEX)  00005320 
C2  =  MIN1F(C2,  YSIZE)  00005330 
IC2=  MAX1F(C2»  SIZEY)  00005340 
C3  =  MIN1F(C3,EXSIZE)  00005350 
IC3=  MAX1FIC3,  SIZEX)  00005360 
C4  =  MIN1FIC4,  YSIZE)  00005370 
IC4=  MAX1FIC4.  SIZEY)  00005380 

244  ICURVU  +  1)  =  IPACK12(  I  CI  .  I  C2  »  I  C3  .  I  C4  )  00005390 

I  I  =  INUM  +  3  00005400 

246  CALL  IPACKLKLABEL.  LABEL1,  IOUMMY)  00005410 
ICURV(II-l)  =  LABEL1  00005420 
ICURV(II)  =  ICURV4  00005430 

9010  IFIMODCURV  -  1)247.247,9015  00005440 

247  CALL  IREADY  (IOUMMY)  00005450 
IF( IDUMMYJ5000, 1260, 5000  00005460 

1260  CALL  IWRITE  (ICONT,  IOUMMY,  1)  00005470 

IF( IDUMMY15C00, 260, 5000  00005480 

260  CALL  IWRITE  (JXTIT,  IDUMMY.ll)  00005490 
IF( IDUMMY)5000, 261 ,5000  00005500 

261  CALL  IWRITE  (JYTIT.  IDUMMY.ll)  00005510 
IF( IDUMMY)5000»265»5000  00005520 

265  IF(  ITD9269. 9268. 9269  00005530 

9268  CALL  IWRITE  (LTITLF.  IDUMMY.  7)  00005540 

IF( IDUMMY)50C0»9269»5000  00005550 

92b9  IF( IT219271, 9270, 9271  00005560 

9270  CALL  IWRITE  (LTITLE(7),  IDUMMY,  7)  00005570 
IF( IDUMMY 15  000,9271,5000  00005  580 

9271  CALL  IWRITE  (KAXIS,  IDUMMY,  5)  00005590 
IF( IDUMMY)5000, 9015, 5000  00005600 

9015  IF( IPOINTJ9020, 270, 9020  00005610 

270  CALL  IWRITE  (ICURV,  IDUMMY,  II)  00005620 

IF( IDUMMY15000, 9020, 5000  00005630 

9020  IFtMODCURV  -  1)272,2  72,^025  00005640 

272  IFIIGRID  -  1)9025.273,9025  00005650 
C          GENERATE  GRID  IF  CALLED  FOR.  00005660 

273  1X100  =  ISIZEX*100  00005670 
IY100  =  ISIZFY»100  00005680 
NEXT1  =  IBOTMGN  00005690 
NFXT2  =  LFTMGN  +  1X100  00005700 
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JGRID(l)  =  0 

DO  1274  J=l,19.2 

JGRID(J+1)  =  IPACK12  (LFTMGN,  NEXT1»  NEXT2.  NEXT1) 

IF(NEXT1  -  IBOTMGN  -  I Y100 ) 127? , 12 75. 1275 

1273 

NEXT1  =  NEXT1  +  100 

JGRID(J+2)  =  IPACK12  (NEXT2.  NEXT1.  LFTMGN.  NEXT1) 

IF(NEXT1  -  IBOTMGN  -  I Y100 ) 1274 . 1276 . 1276 

1274 

NEXT1  =  NEXT1  +  ICO 

1275 

JGRID(J+2)  =  IPACK12  (NEXT2.  M^XTl.  NEXT2 

,  NEXT1) 

1276 

JGRIDU  +  3)  =  ICURV3 
JGRID(J+4)  =  ICURV4 
CALL  IWRITEI JGRID.  I  DUMMY ,  J+4) 
IF( I  DUMMY )5000. 1277,5  000 

1277 

NEXT1  =  LFTMGN 

NEXT2  =  IBOTMGN  +  IY100 

DO  1279  J=1.11.2 

JGRIDIJ+1)  =  IPACK12  (NEXT1,  IBOTMGN.  NEXT1,  NEXT2) 

IFJNEXT1  -  LFTMGN  -  1X100)1278.1280.1280 

1278 

NEXT1  =  NEXT1  +  100 

JGRID(J+2)  =  IPACK12  (NEXT1.  NFXT2.  NEXT1 

♦  IBOTMGN) 

IF(NEXT1  -  LFTMGN  -  1X100)1279,1281.128  1 

1279 

NEXT1  =  NEXT1  +  100 

1280 

JGRID(J+2)  =  IPACK12  (NEXT1,  NFXT2.  NEXT1 

,  NEXT2) 

1281 

JGRID(J+3)  =  ICURV3 

JGRID(J+4)  =  ICURV4 

CALL  IWRITE  (JGRID.  IDUMMY,  J+4) 

IF( IDUMMY) 5 000, 902 5, 50 00 

9025 

IF (I  POINT  19030,276, 90 30 

GENERATE  POINT  PLOT  RECORDS  IF  CALLED 

FOR. 

9030 

IOUT  =  0 

CALL  IPACKL1  (LABEL,  LABEL1,  LABEL2) 

DO  9050  1=1,NUMPTS 

CI  =  X( I )*SCX  +  SHIFTX 

C2  =  Y( I )*SCY  +  SHIFTY 

IFCC1  -  EXSIZE)9031, 9031, 9034 

9031 

IF(C2  -  YSIZEJ9032, 9032, 9034 

9032 

IF(C1  -  SIZFX)9034, 9033, 9033 

9033 

IF(C2  -  SIZEY19034, 9035,9035 

9034 

IOUT  =  IOUT  +1 
GO  TO  9050 

9035 

IC1  =  CI 

GO  TO  (9036.9037.9038.9039.904O) .IPOINT 
GENERATE  CROSS. 

9036 

ICURVI2)  =  IPACK12  (IC1-5,  IC2-5,  IC1+5, 

IC2+5) 

ICURV(3)  =  IPACK12  (IC1   ,  IC2   ,  IC1-5, 

+C2+5) 

ICURVI4)  =  IPACK12  (IC1+5,  IC2-5,  IC1+5, 

IC2-5) 

GO  TO  9041 

GENERATE  PLUS. 

9037 

ICURV(2)  =  IPACK12  (IC1   ,  IC2-5,  IC1   , 

IC2+5) 

ICURVI3)  =  IPACK12  (IC1   ,  IC2   .  IC1-5, 

IC2   ) 

ICURV(4)  =  IPACK12  (IC1+5,  IC2   ,  IC1+5, 

IC2   ) 

GO  TO  9041 

GENERATE  SQUARE. 

9038 

ICURV12)  =  IPACK12  (IC1+4,  IC2-4,  IC1+4, 

IC2+4) 

ICURV(3)  =  IPACK12  (IC1-4,  IC2+4,  IC1-4, 

IC2-4) 

ICURV(4)  =  IPACK12  (IC1+4,  IC2-4,  IC1+4, 

IC2-4) 

00005710 
000C5720 
00005730 
00005740 
00005750 
00005760 
00005770 
00005780 
00005790 
00005800 
00005810 
00005820 
00005830 
00005840 
00005850 
00005860 
00005870 
00005880 
00005890 
00005900 
00005910 
00005920 
00005930 
00005940 
00005950 
00005960 
00005970 
00005980 
00005990 
00006000 
00006010 
00006020 
00006030 
00006040 
00006050 
00006060 
00006070 
00006080 
00006090 
00006100 
00006110 
00006130 
00006140 
00006150 
00006160 
00006170 
00006180 
00006190 
00006200 
00006210 
00006220 
00006230 
00006240 
00006250 
00006260 
00006270 
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GO  TO  9041  00006280 

C          GENERATE  DIAMOND.  00006290 

9039  ICURV12)  =  IPACK12  (IC1+5,  IC2  ,  IC1  ,  IC2+5)  0000630C 
ICURV(3>  =  IPACK12  (IC1-5.  IC2  .  IC1  ,  IC2-5)  00006310 
ICURVJ4)  =  IPACK12  (IC1+5,  IC2  .  IC1+5,  IC2  )  00006320 
GO  TO  9041  00006330 

C          GENERATE  TRIANGLE.  00006340 

9040  ICURVI2)  =  IPACK12  (IC1+5.  IC2-3,  IC1  .  IC2+6)  00006350 
ICURVI3)  =  IPAC<12  (ICl-5,  IC2-3,  IC1+5,  IC2-3)  00U06360 
ICURVI4)  =  ICURV(3)  00006370 

9041  I F ( I  -  NUMPTSJ9043. 9042. 9043  00006380 
904?  ICURV(5>  =  LABFL2  00006390 

GO  TO  9046  00006400 

9043  I F ( I  -  1)9045.9044.9045  00006410 

9044  ICURV(5)  =  LABEL1  00006420 
GO  TO  9046  00006430 

9045  ICURV(5)  =  ICURV3  00006440 

9046  ICURV(6)  =  ICURV4  00006450 
CALL  IWRITE  (ICURV.  I  DUMMY ,  6)  00006460 
IF( IDUMMY15000, 9050, 5000  00006470 

9050  CONTINUE  00006480 

IF< IOUT)9048. 276, 9048  00006490 

9048  PRINT  9104,  IOUT  00006500 

9104  FORMAT  (/,  IX,  12,  29H  POINT(S)  WERE  OFF  THE  GRAPH.  ,/)  00006510 

C          SET  UP  RETURN.  00006520 

276  IF(MODCURV)277,278,277  00006530 

277  IFIMODCURV  -  3)279,278,279  00006540 

278  ITEST  =  0  00006550 
PRINT  130,< JJTITLEI I )  .1  =  1.12)  00006560 

13C  FORMAT  (/,  19H  GRAPH  TITLED  .  .   . 6A8 » / , 19X , 6A8 ,  00006570 

1             24H   .  .  HAS  BEEN  PLOTTED.  ,/,lH0)                       00006580 

IDUMMY  =  ITYP2I IDUMMY)  00006590 

IF( IDUMMY)5670. 656. 5670  00006600 

656  LAST  =  0  00006610 

RETURN  00006620 

279  ITFST  =  1  00006630 
IDUMMY  =  ICLOCM  IDUMMY)  00006640 
LAST  =  0  00006650 
RETURN  00006660 

C          THESE  ARE  THE  NORMAL  RETURNS.  00006670 

C          NOW  SET  UP  THE  RETURN  FOLLOWING  A  TAPE  ERROR.                  00006680 

5000  IFIMODCURV  -  1)5001.5001.5002  00006690 

5001  IDUMMY  =  ITYPEK  IDUMMY  )  00006700 
GO  TO  247  00006710 

5002  PRINT  5100  00006720 
5100  FORMAT  (/.  36H  TAPE  ERROR  IN  WRITING  GRAPH  OUTPUT.  )  00006730 

IDUMMY  =  ITYPEK  IDUMMY  )  00006740 

GO  TO  1007  00006750 

5670  IDUMMY  =  I  TYPE1 ( I  DUMMY )  00006760 

END  00006770 

C  00006780 

C  00006790 

SUBROUTINE  IREADY  'T™'MMYJ  00006800 

C          SELECTS  TAPE  8  (   LL  LOOP  UNTIL  READY).  WRITES  EOF  ON  8.  00006810 

*          MACHINE  LANGUAGE  WILL  NOT  BE  NECESSARY  IN  FORTRAN  62-3.  00006820 

LOCI IFIVE  =5).  00006830 
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FXF  (520413) 

EXF7  (520CCB). 

SELECT  <?EAD  AND  WAIT  TAPE. 

O0006840 

1NEX 

FXF7(0005n3> 

SLJ  (lRfJVJ. 

'»!T  ON  CH  5  ACTIVF. 

00006850 

EXF7(520C"!R) 

SLJ  (1NEX). 

EXIT  ON  TAPE  READY. 

O0006860 

LHA  ( IFIVF ) 

5VJ  (lB'JF). 

T-RMIMATF 

00006870 

13U" 

EXF5(N). 

BUFFER. 

00006880 

1RDY 

EXF  (52041B) 

EXF7( 5200CB). 

SELECT  AND  WAIT  TAPE  8. 

00006890 

ENA  (0). 

CLEAR  A. 

00006900 

EXF  (020003) 

EXF  (520063). 

STOP  CLOCK  AND  BACK-SPACE  8. 

00006910 

EXF7(520013) 

SLJ  (1END). 

EXIT  IF  NOT  AT  LOAD  POINT. 

00006920 

-EXF7(520003). 

WAIT  TAPE  8. 

00006930 

FXF7(520073) 

SLJ  (1E0F). 

EXIT  IF  NO  EOF. 

00006940 

ENA  (I  DUMMY) 

SAU  (2BUF). 

MOVE 

00006950 

INA  (1) 

STA  (IFIVE). 

FORWARD 

00006960 

23UF 

EXF5(N) 

EXF7(5200CB). 

OVER  RECORD. 

00006970 

1E0F 

ENA  (0) 

EXF7(00C613). 

CLEAR  A,  WAIT  CH  6. 

00006980 

EXF  (620413) 

EXF7I620003). 

SELECT  WRITE  AND  WAIT  TAPE. 

00006990 

EXF  (62041B) 

EXF7(62000B). 

SELECT  AND  WAIT  TAPE  8. 

00007000 

FXF  (620033) 

EXF7(620no3). 

WRITE  EOF  AND  WAIT. 

00007010 

EXF7(620073) 

ENA(IO). 

EXIT  ON  NO  END  OF  TAPE. 

00007020 

1END 

STA  (  I  DUMMY). 
END 

00007030 
00007040 

C 

00007050 

SUBROUTINE  I WRI TE (  I  START .  IDUMMY, 

,  IWORDS) 

00007060 

C 

WRITE  RECORD  OF  IWORDS*  STARTING  WITH  ISTART.  PUT  IDUMMY  =  0 

00007070 

C 

IF  RECORD 

CORRECTLY  WRITTEN* 

OTHERWISE  SET  NON-ZERO. 

00007080 

* 

MACHINE  LANGUAGE  WILL  NOT  RE 

NECESSARY  IN  FORTRAN  62-3. 

00007090 

LOC( I  SIX  =  6) . 

00007100 

-FXF7(00061B). 

WAIT  CH  6. 

00007110 

EXF  (620413) 

EXF7(62000B). 

SELECT  WRITEt  WAIT  TAPE. 

00007120 

FXF  (62041P) 

EXF7(620003). 

SFLECT  AND  WAIT  TAPE  8. 

00007130 

ENQ  (1113). 

SET  COUNTER. 

00007140 

1AGN 

ENA  (ISTART) 

I NA ( 1  )  . 

STARTING  ADDRESS. 

00007150 

SAL  (1BUF) 

ADD( IWORDS) . 

TERMINAL  ADDRESS. 

00007160 

1BUF 

STAI IS  IX ) 

EXF6(N). 

SUFFER  OUT. 

00007170 

ENA(O) 

EXF7(62000B). 

CLEAR  A.  WAIT  TAPE  8. 

00007180 

EXF7I62007B) 

SLJ  (1END). 

EXIT  IF  NO  END  OF  TAPE. 

00007190 

FXF7(62003R) 

SLJ  (2AGN). 

FXIT  IF  NO  PARITY  ERROR. 

00007200 

FXF7(62H043) 

SLJ  (2END). 

EXIT  IF  LENGTH  ERROR. 

00007210 

2AGN 

FXF  (620063) 

EXF7(62000B). 

BACKSPACE  AND  WAIT. 

00007220 

QRS  (3) 

QJPK1AGN). 

TRY  WRITE  3  TIMES. 

00007230 

1ENO 

FXF  (620033) 

ENA  (10). 

WRITE  EOF,  NON  ZERO  A. 

00007240 

2END 

STA  (IDUMMY). 
END 

STORE  RESPONSE. 

00007250 
00007260 

C 

00007270 

FUNCTION  ITYP2 

( IDUMMY ) 

00007280 

C 

TYPE  WORD 

GRAPH. 

00007290 

* 

WILL  NEED 

REWRITING  IN  FORTRAN  62-3. 

00007300 

CON(LC  =  57B» 

Ml  =  4513123O15O50CO0B  ). 

00007310 

LOCdTWO  =  2). 

00007320 

-EXF7I00061B) . 

WAIT  CH  6. 

00007330 

FXF  (62041R) 

EXF7(620003). 

SELECT  AND  WAIT  TAPF. 

00007340 

FXF  (62041B) 

FXF7(62nOOB). 

SFLECT  AND  WAIT  TAPE  8. 

00007350 

FXF  (620033) 

EXF  (010003). 

WRITE  EOF.  START  CLOCK. 

00007360 

ENA  (0) 

EXF7(62000B). 

CLEAR  A.  WAIT  TAPE  8. 

00007370 

FXF7I62007B) 

ENA  (10). 

EXIT  IF  NO  END  OF  TAPE. 

00007380 

STA  (ITYP2). 

STORE  RESPONSE. 

00007390 
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-EXF7(00021B).  WAIT  CH  2.                        00007400 

EXF7(11141B)     SLJ  (1TYP).  EXIT  IF  UPPER  CASE.               00007410 

EXF  (211O0B)     ENA  (LC+1).  TYPE                              00007420 

STA  (ITWO)       EXF2(LC).  LOWER  CASE.                    00007430 

-EXF7(0C021B).  WAIT  CH  2.                        00007440 

1TYD   FXF  (211008)     ENA  (Ml+1).  TYPE                              00007450 

STA  (ITWO)       EXF2(M1).  GRAPH                           00007460 

END  00007470 

C  00007480 

FUNCTION  ITYPE1  (I  DUMMY)  00007490 

C          REWIND  TAPE  8.  REQUEST  NEW  TAPE.  AND  WAIT  TILL  READY.  00007500 

*  WILL  NEED  REWRITING  IN  FORTRAN  62-3.  00007510 
CON(LC  =  57B»  Ml  =  4 5 1 5 1 12030242004B  .  M2  =  1 10330220406203 1 B ,  00007520 

1     M3  =  0401301520043342B).  00007530 

RSVtMESS  =  3).  00007540 

LOCUTWO  =  2).  00007550 

-EXF7I00061B) .  WAIT  CH  6.                        00007560 

EXF  (62041B)     EXF7 ( 62000B ) .  SELECT  AND  WAIT  TAPE.            00007570 

EXF  (62041B)     EXF7 ( 62000B ) .  SELECT  AND  WAIT  TAPE  8.          00007580 

EXF  (62003B)     EXF7  { 62000(3 )  .  WRITE  EOF  AND  WAIT.               00007590 

EXF  (62007B)     EXF7 ( 0002  IB ) .  REWIND  WITH  INTERLOCK,  WAIT  CH  200007600 

EXF7(11141B)     SLJ  (1TYP).  EXIT  IF  UPPER  CASE.               00007610 

EXF  (21100B)     ENA  (LC+1).  TYPE                              00007620 

STA  (ITWO)       EXF2ILC).  LOWER  CASE.                     00007630 

1TYP   LDA  (Ml)          STA  (MESS).  00007640 

LDA  (M2)         STA  (MESS+1).  00007650 

LDA  (M3>         STA  (MESS+2).  00007660 

-EXF7I00021B).  WAIT  CH  2.                        00007670 

EXF  (21100B)     ENA  (MESS+3).  TYPE                              00007680 

STA  (ITWO)       EXF2(MESS).  MESSAGE.                        00007690 

-EXF7(00061B).  WAIT  CH  6.                       00007700 

-EXF7(62000B) .  WAIT  TAPE.                        00007710 

EXF  (62041B)     EXF7 ( 62000B ) .  SELECT  AND  WAIT  TAPE  8.          00007720 

EXF  (010008).  START  CLOCK.                     00007730 

END  00007740 

C  00007750 

FUNCTION  ICODE  (ANUMBER)  00007760 

C          CODES  ABSOLUTE  VALUE  OF  A  FLOATING  POINT  NUMBER  (BETWEEN  00007770 

C          1.0E-100  AND  1.0E+100)  INTO  8-CHARACTER  BCD  WORD  OF  THE  FORM   00007780 

C          1.23E+45.  ICODE  =  8H0.00E+00  IF  MAGNITUDE  OUT  OF  RANGE.  00007790 

*  CHECK  AVAILABILITY  OF  LIBRARY  FUNCTIONS  IN  FORTRAN  62-3.  00007800 
DIMENSION  11(8)  00007810 
BNUMBER  =  ABSF( ANUMBER)  00007820 
IF(BNUMBER  -  1 . OE+100 ) 7 .6 » 6  00007830 

7  IFIBNUMBER  -  1 . 0E- 100 ) 6 »6 » 2  00007840 

6  ICODE  =  8H0.00E+00  00007850 

RETURN  00007860 

C          THIS  IS  ERROR  EXIT.  00007870 

2  CALL  SCALEIT  (BNUMBER.  ISCALlO.  FACTOR.  3)  00007880 
ISIGNSC  =  XSIGNFf 1. ISCALlO)  00007890 
ISCALlO  =  XABSFt ISCALlO)  00007900 
IFACT  =  FACTORMOO.001  00007910 
11(8)  =  XMODF( ISCALlO, 10)  00007920 
II(7)=ISCAL10/10  00007930 
IF(ISIGNSC)4.3,3  00007940 

3  II  (6)  =  8H        +  00007950 
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G  J  TO  5 

I  I  (6)  =  8H 

I  I  (5)  =  8H        E 

I  I  (4)  =  XMODFt IFACT.iO ) 

I  I  (3)  =  fXMODFt !FACT»100) J/10 

I  I  (2  )  =  8H        . 

11(1)  =  IFACT/IOO 

CALL  IPACK  (II,  IPACKFD) 


ICODc  = 

RFTURN 

FND 


IPACKFD 


SUBROUTINE  SCALEIT  ( ANU^BER . I SCAL 10 ,F AC TOR .MODE > 
C  FINOS  FACTOR  (BETWEEN  1.0  AMD  9.99...)  AND  SCALE  OF 

C  DEFINED  BY      ANUMBER  =  FACTOR* 1 C. ** I SCAlIO. 

C  MODE  IS  THE  NUMBER  OF  SIGNIFICANT  FIGURES  REQUIRED. 

C  BE  BETWEEN  1  AND  10  OR  IT  WILL  BE  DUT  EQUAL  TO  SIX. 

*  CHECK  AVAILABILITY  OF  LOGl^F  IN  FORTRAN  62-3. 

ISCAL10=LOG10F(ANUV.BFR) 

FACTOR  =  ANUMBER/10.**ISCAL1C- 

IFfFACTOR  -  t.l)l,2,2 

1  FACTOR  =  FACTOR*100. 
1SCAL10  =  ISCAL10  -  2 
GO  TO  8 

2  IFIFACTOR  -  1.0)3.8.4 

3  FACTOR  =  FACT0R*10. 
ISCAL10  =  ISCAL10  -  1 
GO  TO  8 

4  IFfFACTOR  -  100.0)6.5.5 

5  factor  =  factor/ioo. 
iscal10  =  iscal10  +  2 

GO  TO  8 

IF(FACTOR  -  10.0)8*7,7 

FACTOR  =  FACTOR/10. 

ISCAL10  =  ISCAL10  +  1 

IF (M0DE)9»9,10 

MODE  =  6 

GO  TO  11 

IFJMODE  -  10)11,11,9 

IFACTOR  =  FACT0R*10.**(M0DE  -  1)  +  0.5 

FACTOR  =  IFACTOR 

FACTOR  =  FACTOR/10. **<MODE  -  1) 

IF(FACTOR  -  10.  J13.12.12 

FACTOR  =  1. 

ISCAL10  =  ISCAL10  +  1 

RFTURN 

END 


6 

7 

8 
9 

10 
11 


12 


13 


SUBROUTINE  [PACK  (II.  IPACKED) 
C  TAKES  8  SIX-BIT  WORDS  AND  PACKS  THEM  LEFT  TO  RIGHT 

C  IN  IPACKED.  IF  WORD  IS  ZERO.  12B  IS  SUBSTITUTED. 

*  CONVERT  TO  CODAP  FOR  FORTRAN  62-3. 

CON( IZERO  =  12B). 
SIUKI5AVE)      ENIK8). 
1NFX   LDAKII)         AJP1(2NEX>. 
LDA  ( IZERO). 


" T007560 
03C07970 
0C0C798C 
C0C37990 
"0008000 
C0008010 
C000802C 
0C0C8O3O 
00008040 
00008050 
00008060 
00008070 
00008080 

10  AS  00008090 
00008100 

THIS  MUST  00008110 
00008120 
00008130 
00008140 
00008150 
00008160 
00008170 
00008180 
00008190 
00008200 
00008210 
00008220 
00008230 
00008240 
00008250 
00008260 
00003270 
00008280 
00008290 
00008300 
00008310 
00008320 
00008330 
00008340 
00008350 
00008360 
00008370 
00008380 
00008390 
00008400 
00008410 
00008420 
00008430 
00008440 
00008450 
00008460 
00008470 
00008480 
00008490 
00008500 
00008510 
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2NEX 


1NFX 


2NEX 


LRS  (6) 
ISKK-1) 
STQ  ( IPACKED) 
END 


SUBROU 

IN 

WIL 

WA 

CON( 13 
SIUK 
ENA  ( 
LOOK 
STAK 
ISKK 
LDO  ( 
ENIK 
ENA  ( 
LDQK 
INIK 
INIK 
ISKK 
LDQ  ( 
ENIK 
LIUK 

END 


TINE  ISHI 

SERTS  026 

L  HAVE  TO 

TCH  ARRAY 

LANK  =  20 

ISAVE  ) 

2). 

ITITLE) 

LTITLE) 

6) 

IBLANK) 

7) 

2). 

ITITLE) 

1) 

-1) 

12) 

IBLANK) 

14) 

ISAVE). 


INI  K-2). 
SLJ  (1NEX). 
LIU1 ( ISAVE). 


FT6  (  ITITLE,  LTITLE) 

AHEAD  OF  6-WORD  TITLE  RECORD. 

BE  CONVERTED  TO  CODAP  IN  FORTRAN  62-3. 

INDEXING  IN  FORTRAN  62-3. 
202O202O2O2C20R). 
EMI  1(1). 


SET  COUNT, 


LLS 
LLS 
SLJ 
LLS 


(42)  . 
(6). 

( 1NEX). 
(42)  . 


STAKLTITLE). 

LLS  (42). 
STAKLTITLE). 
LLS  (6). 
SLJ  (2NEX). 
LLS  (42). 
STAKLTITLE). 


SAVE  INDEX. 
ENTER  02B. 
PERFORM 

SHIFTING. 
CHECK  IF  COMPLETE. 
COMPLETE  LAST  WORD. 
STORE  LAST. 
REPEAT 
FOR 

SECOND 
TITLE 
LINE. 


RESTORE  INDEX. 


FUNCTION  IPACK12  (  I  ONE • I  2 » I  3 . I  4 ) 

PACKS  FOUR  12-BIT  WORDS  INTO  ONE  48-BIT  WORD. 
WILL  REQUIRE  CONVERSION  TO  CODAP  IN  FORTRAN  62-3, 


LDA( IONE) 

QLS(36) 

LDQ( 13) 

LLS(12) 

QLS(36) 

STA(  IPACK12  ). 

END 


LDQ( 12) . 
LLS(12) . 
OLS(36) . 
LDQ( 14) . 
LLS( 12)  . 


SUBROUTINE  IPACKL1  (LABELt  LABFL1,  LABEL2) 
PACKS  TWO  4-CHARACTER  LABFLS. 
USE  DECODE/ENCODE  IN  FORTRAN  62-3. 

CON( IFLAG  =  37773777B)  . 


LLS(24) 
LDA  ( IFLAG) 
LDA  (  IFLAG) 
STA  (LABEL2). 
END 


STA  (LABEL1). 
LDQ  (LABEL). 
LLS(24). 


FUNCTION  ICLOCKI IDUMMY) 

FXF  (01000B). 
END 

IC2  =  C2 
END 
END 


START  CLOCK, 


C0008520 
00008530 
00008540 
00008550 
00008560 
00008570 
00008580 
00008590 
00008600 
00008610 
00008620 
00008630 
00008640 
00008650 
00008660 
00008670 
00008680 
00008690 
00008700 
00008710 
00008720 
00008730 
00008740 
00008750 
00008760 
00008770 
00008780 
00008790 
00008800 
00008810 
00008820 
00008830 
00008840 
00008850 
00008860 
00008870 
00008880 
00008890 
00008900 
00008910 
00008920 
00008930 
00008950 
00008940 
00008960 
00008970 
00008980 
00008990 
00009000 
00009010 
00009020 
00006120 
00009030 
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APPENDIX  B 
PROGRAM  NOMENCLATURE 
The  following  listing  of  program  symbols  is  provided  to  help  the 
reader  understand  symbols  used  in  the  digital  computer  program.  Symbols 
identified  in  section  3  as  well  as  symbols  readily  identifiable  are  not 
listed  in  this  appendix. 


X(Y) 

T(D 

AZ 

ANG(L) 

SMAX(L) 

SMIN(L) 

THETA(L) 

SIGXX(LL) 

TRUXY(L) 

PHI(L) 

XSXX(L) 
YSXX(L) 

XTXY(L) 
XTXY(L) 

zx(i,i) 

zx(7,D 
11(1,1) 

XY(7,1) 

zzz 

PHI  1(1) 
SIGN  7(1) 


individual  strain  gage  readings 

percentage  of  depth  of  individual  strain  gage  readings 

the  zero  depth  intercept  of  strain  axis 

the  angle  of  (Ym*  from  gage  leg  number  1 

GVun 

eP 

experimental  value  of  (Tn 

experimental  value  of  T*y 

the  angle  0 

1  components  of  experimental  values  of  <Tr\  in  proper  form 
J  for  plotting 

1    components  of  experimental  values  of  Xxy  in  proper  form 
J  for  plotting 

>  seven  coefficients  used  in  computing  smoothed  value  of  On 
\    seven  coefficients  used  in  computing  smoothed  value  of  txy 


the  angle  in  radians  on  the  circumference  of  the  piping  between 
two  consecutive  smoothed  values  of  On  or  T*y 

the  angle  in  radians  from  rosette  number  12  of  a  smoothed 
value  of  Oiig  ortxY 

smoothed  value  of  On . 
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XSNR(I 
YSNR(I 


TAUN  7(1)       smoothed  value  of  Tky 

V    components  for  plotting  zero  stress  circle  on  On  plot 

I    components  for  plotting  zero  stress  circle  on  T*y  plot 
1  components  of  smoothed  value  of  Gm  in  proper  form  for 


XTNR(I 
YTNR(I 

XSN7(I 
YSN7(I 

XTN7(I 
YTN7(I 


plotting 


~)   components  of  smoothed  value  of  txy  in  proper  form  for 
J  plotting 


72 


APPENDIX  C 

DESCRIPTION  OF  SUBROUTINE  GAUSS  2 

A.  IDENTIFI CATION 

TITLE:  Solution  of  Simultaneous  Linear  Algebraic  Equations 

CO-OP  ID:  F2 

FORTRAN:  Gauss  2 

CATEGORY:  Simultaneous  Linear  Equations 

PROGRAMMER:  C.  B.  Bailey  and  Mary  Haynes 

DATE  OF  COMPLETION:  April  1963 

B#  PURPOSE:  To  solve  one  or  more  sets  of  linear  algebraic  equations  using 
Gaussian  elimination  with  row  pivoting  and  back  substitution, 

C.  USAGE: 

1.  Calling  sequence 

Call  Gauss  2  (N,  M,  EP,  A,  B,  X,  KER) 

N  -  order  of  the  matrix  in  the  equation  Ax  =  b,  that  is,  the 

number  of  linear  equations  (maximum  =  50)  • 

M  -  number  of  vectors  b  for  which  solution  vectors  x  are  to  be 

obtained,  that  is,  the  number  of  sets  of  linear  equations  (Maximum  =  60)« 

EP  -  matrix  condition  parameter  (see  Mathematical  Method). 

A  -  the  elements  of  the  matrix  of  coefficients  of  the  equations, 

stored  in  the  form  (A^i) 

i  =  1,  2,  ...,  N 

J  =1,  2,  ...,  N 

A  is  dimensional  as  a  50  x  50  array. 

B  -  the  elements  of  the  column  vectors  b  etc.,  stored  in  the 

form  (B3J 

i  =  1,  2,  ...,  N 

3=1,  2,  ...,  M 

B  is  dimensional  as  a  50  x  60  array. 

X  -  the  components  of  the  x  vectors  stored  in  the  form  (Xjj) 
i  -  1,  2,  ...,  N 

j  =  1,  2,  ...,M 
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X  is  dimensional  as  a  50  x  60  array, 

KER  -  error  flag  (See  Error  Return) 
2.  Space  Required: 
3«  Temporary  or  Common  Storage  Required 

4.  Error  Return:  Argument  KER 

1.  indicates  no  errors 

2,  indicates  that  matrix  is  singular  or  nearly  singular. 

5,  Accuracy:  Not  applicable. 
D.  METHOD 

This  subroutine  is  an  adaptation  of  the  program  F2  VTEX  LINEQN. 
Please  consult  the  write-up  of  that  program  for  the  mathematical  method, 
The  subroutine  gives  no  printout. 
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APPENDIX  D 
DESCRIPTION  OF  SUBROUTINE  DRAW 

A.  IDENTIFICATION: 

TITLE:     General  Graph  Output  Subroutine 

CO-OP  ID:     J7-NPS-DRAW 

CATEGORY:     Output  for  Off -Line  Plotting 

PROGRAMMER:     J,  R.  Ward 

DATE:     February  1964}  REVISED  JUNE  1965 

B.  PURPOSE: 

This  subroutine,  when  provided  with  the  necessary  information,  generates 
a  magnetic  tape  in  the  proper  format  for  subsequent  off-line  graph  plotting 
using  the  CDC  160  or  160A  GRAFPLOT  program  and  a  CalComp  165  Plotter 
(see  references  1  and  2).  Provision  is  made  for  curve  drawing  and  point 
plotting,  automatic  scaling,  graph  titling  and  axis  annotating.  An 
attempt  was  made  to  provide  a  considerable  amount  of  flexibility,  at  the 
expense,  necessarily,  of  a  relatively  large  number  of  arguments  and  a 
rather  high  memory  requirement. 

C.  USAGE: 

1.  Definitions: 

In  what  follows  the  word  "graph"  will  be  taken  to  mean  one  piece 
or  frame  of  graph  paper  on  which  there  may  be  plotted  one  or  more  curves 
and/or  sets  of  points.  A  "curve"  will  mean  a.  continuous  line  generated 
by  the  sequence  of  straight  lines  joining  successive  points  of  the  set 
defining  the  curve.  A  "point  plot"  will  describe  the  representation  of 
a  succession  of  points  by  means  of  symbols  (such  as  a  cross)  on  the  graph. 
The  points  are  not  connected  in  a  point  plot. 

2.  Calling  Arguments: 

All  necessary  information  is  transferred  to  DRAW  through  the  calling 
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arguments.     The  call  statement  is:     CALL  DRAW  ( NUMPTS, X,Y,M0DCURV, 

ITYPS,  LABEL,  ITITLK,  EXSCALK,  YSGALE,  IXUP,  IYRIGHT,MODEXAX,MODhTAX,  IWIDE, 

IHIGH,IGRID,LAST) 

It  is  important  to  realize  that  one  and  only  one  curve  or  set  of  points 

is  plotted  each  time  DRAW  is  called.     However,  it  is  possible  to  call 

DRAW  repeatedly  if  several  curves  and/or  sets  of  points  are  wanted  on 

one  graph. 

The  calling  arguments  are  as  follows: 

a.  NUMPTS:  The  number  of  points  defining  a  curve  (2  ^  NUMPTS 

<  900),  or  the  number  of  points  to  be  point  plotted 
(2   <    NUMPTS  <    30). 

b.  X:  The  array  of  X-ordinates  (  JX^I    ^  10°°  for  i=l,  2,..., 

NUMPTS).     X  must  be  dimensioned  at  least  equal  to 
NUMPTS  in  the  calling  program. 

All  points  will  be  considered  to  have  the  same 
X-ordinate  if      |  X^^    -    ^min'    ^*  10""°'.     The  common 
value  will  be  put  equal  to  zero  if       |  Xjnax  I   ^  lO-0?. 

c.  Y:  The  array  of  Y-ordinates,  with  properties  corresponding 

to  the»X-ordinates,  above.     Y  must  be  dimensioned  at 

least  equal  to  NUMPTS  in  the  calling  program. 

Controls  the  number  of  curves,  and/or  sets  of  points 

on  one  graph: 
=  0     •  This  is  the  only  curve,  or  set  of  points,  to  be  plotted 

on  this  graph. 
=  1  This  is  the  first  of  two  or  more  curves,  and/or  sets  of 

points,  to  be  plotted  on  this  graph. 
=  2  This  is  an  intermediate  curve,  or  set  of  points. 

=  3  This  is  the  last  curve,  or  set  of  points,  for  this  graph. 


d.     M0DCURV: 
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e.  ITYPE:       Controls  the  type  of  plot  (i.e.,  curve  or  point  plot): 

■  0     This  set  of  points  is  to  be  represented  by  a  curve. 
=  1    These  points  are  to  be  plotted  with  a  cross  (x). 

■  2    These  points  are  to  be  plotted  with  a  plus  (+). 

■  3     These  points  are  to  be  plotted  with  a  square  (  □  ). 
=  4    These  points  are  to  be  plotted  with  a  diamond  (  O  ). 

■  5     These  points  are  to  be  plotted  with  a  triangle  (  ^  ). 

f .  LABEL:       This  is  a  Hollerith  curve  or  point  identifier.  If  a 

curve  is  being  drawn,  LABEL  must  have  4  characters 
(including  any  blanks),  and  these  will  be  reproduced 
beside  the  end  of  the  curve.  This  argument  can  be  set 
in  the  calling  program  by  a  statement  such  as 

LABEL  =  4H  ONE         (A=  blank) 
or   LABEL  =  4H1234 
or   LABEL  =  4HAAAA      The  latter  must  be  used  when 

^  AAAA  ■   i 

no  label  is  wanted. 

If  a  set  of  points  is  being  plotted,  LABEL  is  an  8- 
character  identifier.  The  first  4  characters  will  be 
reproduced  beside  the  first  point,  and  the  last  four 
characters  will  be  reproduced  alongside  the  last  point. 
This  argument  can  be  set  by  statements  such  as 

LABEL  =  8HFRSTLAST 
or    LABEL  =  8R\,AAA  ONE 

AAAA 

or    LABEL  =  8HA0NEA123 


or 


LABEL  =  8HAAAAAAAA     The  latter  must  be  used  if 

A/nAAAAAA  i— ■■ 


no  labels  are  wanted. 

The  above  arguments,  a.  through  f.  (and  q.),  have 

meaning  every  time  DRAW  is  called.  On  the  other  hand, 
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g.     ITITLE: 


h.     EXSCALE: 


f.     YSGALE: 


j.     IXUP: 


k.     IYRIGHT: 


1.     MODEXAX: 
=  0 


the  remaining  arguments,  g.  through  p.,  have  no 
meaning  except  when  MODCURV  =  0  or  1. 
An  array  of  twelve  8-character  Hollerith  words,  the 
first  six  of  which  will  form  the  first  title  line,  and 
the  last  six  the  second.     The  array  must  be  dimensioned 
12  in  the  calling  program,  must  contain  the  user's  job 
identification,  and  must  have  unwanted  characters  set 
to  blank.     For  example: 

DO  1  I  =  1,12 
1      ITITLB(I)  =  8H 

ITITLE(l)  =  8HASMITH,A 

ITITLE(2)  -  8HJ.AJ.AAA 

ITITLE(7)  =  8HATESTIT. 
X-scale  in  units  per  inch  (10""  <  EXSCALE  <   1099). 
EXSCALE  will  always  be  rounded  off  to  one  figure  sig- 
nificance.    If  EXSCALE  =  0,  the  X-scale  will  be  computed 
by  DRAW.     This  is  called  auto-scale. 

Y-scale  in  units  per  inch,  with  properties  corresponing 
to  those  of  EXSCALE. 

Distance,  in  inches,  of  the  X-axis  from  the  bottom  of 
the  graph  (0  <  IXUP  <  IHTGH).       This  argument  will  be 
ignored  unless  MODEXAX  =  2,   see  below. 
Distance,  in  inches,  of  the  Y-axis  from  the  left  of  the 
graph  (0  ^  IYRIGHT  ^  IWIDE).     This  will  be  ignored  unless 
MODEYAX  =2,   see  below. 

Determines  the  mode  of  the  X-axis  location: 
The  X-axis  will  be  located  automatically  by  DRAW,  with 
the  origin  of  Y  on  the  graph. 
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m.  MODEYAX: 


n.  IWIDE: 


o.  IHIGH: 


IGRID: 


=  1      The  X-axis  will  be  automatically  located  by  DRAW, 

with  the  origin  of  Y  removed  (in  one's  imagination)  an 
integer  number  of  inches  above  or  below  the  graph,  if 
this  is  appropriate.  This  option  can  be  used  only  if 
the  Y-scaling  is  automatic  (YSCALE  =  0). 
=  2      The  X-axis  location  will  be  as  specified  by  IXUP. 

determines  the  mode  of  Y-axis. location  in  the  same  way 
as  MODEXAX,  above,  governs  the  X-axis  location. 
Width  of  graph  in  inches  (l  <   IWIDE  £  9).  If  IWIDE 
is  out  of  this  range,  a  value  of  8  will  be  assumed. 
Height  of  graph  in  inches  (l  ^  IHIGH  ^15).  If  IHIGH 
is  out  of  this  range,  a  value  of  8  will  be  assumed. 
If  IGRID  =  1,  a  1'*  x  1"  grid  will  be  superimposed  on 
the  graph.  This  is  useful  only  if  plain  paper  is  used 
on  the  CalComp  Plotter. 

Indicates  to  the  calling  program  whether  the  previous 
plot  was  completed  successfully.  The  codes  are: 
Last  plot  was  completed  successfully. 
Last  plot  was  not  completed  successfully. 
Last  plot  was  not  completed  successfully,  and  no  further 
graph  output  will  be  attempted  until  DRAW  is  next 
•  entered  with  MODCURV  =  1  or  0. 

An  attempt  was  made  to  enter  DRAW  with  MODCURV  /  1  or 
0  while  the  error  lockout  was  set. 

This  argument  must  always  be  a  variable  name  and  never 
a  number  in  the  call  statement. 
3.  Notes  and  Comments? 

a.  The  graph  scales-and,  if  MODEXAX  =  1  and/or  MODEYAX  =  1,  the 


LAST: 


=  0 

=  1 
=  2 


=  3 
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amounts  of  origin  offset  -  are  always  output  as  part  of  the  graph  title. 

b.  Each  time  a  graph  is  completed,  a  message  to  this  effect  is 
printed  on  both  the  standard  output  and  the  console  typewriter. 

c#  There  are  internal  checks  of  the  input  to  DRAW  to  prevent  incorrect 
use.  If  an  input  error  is  detected,  an  attempt  will  be  made,  where 
possible,  to  complete  the  plot.  If  an  argument  is  "corrected"  in  this 
process,  the  user  will  be  so  informed  on  the  standard  output.  If  it  is 
not  possible  to  complete  the  plot,  the  user  will  be  informed  of  the  reason 
by  a  message  on  the  standard  output. 

d.  If  part  or  all  of  a  curve  would  fall  more  than  0.6"  laterally 
beyond  the  ends  of  the  X-axis,  or  0.7"  vertically  beyond  the  ends  of  -the 
Y-axis,  the  X  and/or  Y  ordinates  will  be  limited  so  that  the  curve  will 

.  typically  become  a  line  along  part  or  all  of  the  boundary  of  the  graph 

as  here  defined. 
I 

e.  If  one  or  more  points  of  a  point  plot  would  fall  outside  the  graph 
area,  the  plot  of  that  point,  or  poirfts  will  be  inhibited.  The  number  of 
such  points  will  be  reported  to  the  user  on  the  standard  output. 

f .  It  should  be  pointed  out  that  the  X  and  Y  scaling  and  axis  locating 
processes  are  entirely  independent,  so  that,  for  example,  X  might  be 
auto-scaled,  while  the  Y-scale  is  specified.  At  the  same  time  the  X-axis 
might  be  located  automatically,  while  the  Y-axis  location  is  specified. 

g.  It  must  be  remembered  that  the  scales  and  axis  locations  of  a 
multi-plot  graph  are  set  when  DRAW  is  called  for  the  first  time  (with 
MODGURV  =1).  Thus  the  user  must  attempt,  at  that  time,  to  achieve 
scaling  and  axis  location  which  will  be  appropriate  to  all  the  plots  he 
intends  to  make  on  the  one  graph.  Particularly  if  the  automatic  features 
of  DRAW  are  selected,  foresight  will  be  demanded  of  the  user  in  this 
respect. 
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4*  Auto-Scale  Properties: 

The  scale  factor  is  chosen  from  amongst  the  values  1,2,  or  5  units 
per  inch,  or  some  power  of  10  times  one  of  those  values.  A  curve,  or 
set  of  points  which  is  plotted  with  auto-scale  will  normally  lie  entirely 
within  the  graph  area  as  defined  in  3.d.,  above.  The  only  exception  may 
occur  if  an  axis  is  placed,  by  the  user,  along  one  edge  of  the  graph 
(e.g.,  IXUP  =  0,  MODEXAX  =  2).  In  such  a  case,  points  "outside"  the  axis 
are  not  considered  in  the  selection  of  a  scale  factor  (e.g.,  negative  Xj[ 
do  not  affect  the  choice  of  scale  when  IXUP  =0).  If  automatic  axis 
location  as  well  as  auto-scale  is  selected,  the  plot,  if  it  does  not  fill 
the  graph  area,  will  be  placed  as  far  as  possible  towards  the  bottom-left 
of  the  graph  area  consistant  with  the  fact  that  the  axes  can  be  set  in 
increments  of  1"  only. 

5.  Space  Required:  3960  cells  (excluding  the  input  arrays). 

6.  Temporary  Storage:  None. 

7.  Error  Print-Outs:  There  are  a  large  number  of  possible  error  print- 
outs. These  are  all  self-explanatory. 

8.  Error  Returns:  All  error  returns  are  preceded  by  self-explanatory 
error  printouts.  An  error  indication  is  transferred  back  to  the  calling 
program  through  the  argument  LAST. 

9.  Error  Stops:  None. 

10.  Tape  Mountings:  Logical  Tape  #8  will  receive  the  binary  graph  output. 
The  standard  monitor  output  will  receive  the  messages  to  the  user. 

11.  Output  Format:  The  format  of  the  binary  graph  output  records  on 
magnetic  tape  is  described  in  reference  1.  The  only  difference  is  that 
in  this  program  the  interpolation  option  is  by-passed  (Set  to  zero  in 
the  graph  output  record).  See  reference  2. 

12.  Selective  Jump  and  Stop  Settings:  None. 
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13.  Tinning:  Variable,  depending  upon  the  number  of  points  and  the 
options  chosen.   Typically  less  than  one  second  per  curve  or  point  plot, 
14#  Accuracy:  The  accuracy  of  results  is  equal  to  the  resolution  of 
the  CalComp  Plotter,  that  is,  0.01"  in  both  the  X  and  Y  directions. 
15.  Equipment  Configuration:  CDC  1604  with  FORTRAN  60  compiler  and 
Library.  A  CDC  160  or  160A  with  CalComp  165  Plotter  is  needed  for  the 
off-line  plotting  using  the  appropriate  GRAFPLOT  program. 
D.  REFERENCES: 

1.  Weir,  Maurice  D.,  Spritzer,  Milton  and  Mcllhenny,  D.W.,  rtl60-A 

Graph  Plot  Program,  ■  Ident  *B001,  SWAP  Library,  15  August  1962. 

2.  Hogg,  R.  L.  and  Glover,  D.C.,  "160  Graf plot  Routine",  Writeup 

Available  from  Computer  Facility,  U.  S.  Naval  Postgraduate 
School,  1  April  1963. 

3.  U.S.  Naval  Postgraduate  School,  Thesis,  "Control  System 

Programming,  Remote  Computing  and  Data  Display,"  by  Robert  Lee 
Hogg  and  Dennis  C.  Glover,  1963. 

« 

N,B.  References  1  and  2  are  included  in  Reference  3  as 
Enclosures  2  and  1,  respectively. 


82 


APPENDIX  E 
PLOTS  OF  DOCK  TRIAL  DATA 
Scale  of  stress  has  been  omitted  for  reasons  of  security 
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FIGURE  E-l 

Polar  plot  of  On  for  rosette  ip.ng  #1 
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FIGURE  E-2 
Polar  plot  of  txy  for  rosette  ring  #1 
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FIGURE  E-3 
Polar  plot  of  On  for  rosette  ring  #2 
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FIGURE  E-4 

Polar  plot  of  t*vf  for  rosette  ring  #2 
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FIGURE  E-5 
Polar  plot  of  On  for  rosette  ring  #3 
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FIGURE  E-6 
Jftftfrr  plot  *f  tkY  for  rosette  ring  #3" 
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FIGURE  E-7 
Polar  plot  of  On  for  rosette  ring  #4 
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FIGURE  E-8 

Polar  plot  of  %cyfor  rosette  ring  #4 
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FIGURE  E-9 
*  Pblar  plot  of  CTn  for  rosette  ring  #5 
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FIGURE  E-10 
Polar  .plot  of  t*Y  for  rosette  ring  #5 
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FIGURE  E-ll 


Polar  plot  of  GTs<  for  rosette  ring  #6 
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FIGURE  E-12 
Polar  plot  of  txy  for  rosette  ring  #6 
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FIGURE  E-13 
Polar  plot  of  (T*  for  rosette  ring  #7 
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FIGURE  E-14 
Polar  plot  of  txY  for  rosette  ring  #7 
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FIGURE  £-15 
Polar  plot  of  QVj  for  rosette  ring  #8 
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FIGURE  E-16 

Polar  plot  of  ^y  for  rosette  ring  #8 
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